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1. Introduction. It is not generally recognized that physical re- 
search is quite as dependent upon aids to the mind as it is upon 
aids to the senses. Indeed, it is the aids to the mind that are fashioned 
with greatest difficulty; only at long intervals does there appear a 
figure possessing that high degree of originality necessary for the con- 
struction of a new mental device that may be used successfully in 
analyzing the complex phenomena of nature. In chemistry, we have 
three such devices, namely: thermodynamics, statistical mechanics, 
and quantum mechanics. For the first two of these we are largely 
indebted to Josiah Willard Gibbs. Although a mathematician by 
training, Gibbs possessed that intuitive physical sense which enabled 
him to reduce physical phenomena to their lowest terms. His con- 
tributions to chemical science rank in importance with those of La- 
voisier and Dalton. 

To the layman it must appear that the rapid progress of modern 
science is due to experimentation. Nothing could be farther from the 
truth. In all ages, there have been men who experimented, just as 
there have been those who philosophized; but their contributions to 
scientific knowledge were few in number and limited in scope. What 
distinguishes our present age from all others is that experiment and 
theory have been joined in complementary fashion. The achieve- 
ments of science have been due to this happy combination. 

The subject that I propose to discuss is a small segment of a much 
broader one, namely, the relation between electricity and matter. 
While the problem of electrolytes is but a fraction of the broader 
problem, its study has led to the development of concepts which have 
made it possible to interpret matter in all its forms in terms of elec- 
trically charged particles. f 

Interesting and instructive as it might be, an historical introduc- 
tion is not possible on this occasion. Suffice it to say that the problem 


* The fourteenth Josiah Willard Gibbs Lecture, delivered at Indianapolis, Indi- 
ana, December 28, 1937, under the auspices of the American Mathematical Society, 
at a joint meeting of the Society with the American Association for the Advancement 
of Science. 

+ This excepts the neutron, the relation of which to other types of ultimate par- 
ticles of matter is as yet little understood. 
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of electrolytes had its inception in the observations of chemists, sev- 
eral hundred years ago, and that on the physical side it received its 
first impetus through the discoveries of Volta. These were followed 
by the epoch-making investigations of Faraday, who laid the founda- 
tions of the theory of electrolytes. He showed that matter and elec- 
tricity are related in a definite, quantitative manner, whatever the 
nature of the atoms concerned. Most important, he introduced the 
concept of an ion as a charged atom or molecule. He accounted for 
the passage of electricity through a solution of an electrolyte by the 
motion of oppositely charged ions through the solution. 

The ionic concept is one of the most important that has ever been 
introduced into physical science. It is more fundamental than the 
concept of atoms and molecules, for we now know that all atoms are 
composed of positively charged nuclei, surrounded by negative elec- 
trons equal in number to the charge on the nucleus. In final analysis, 
all matter is constructed of assemblies of ions. Today, the structure 
of the nucleus itself is interpreted in terms of ionic phenomena. 

The application of the ionic concept to the problem of the structure 
of matter is a by-product of the problem of electrolytes. This latter 
problem is still far from solution in complete form. Through the work 
of Clausius, Hittorf, Kohlrausch, Arrhenius, and others, it finally be- 
came clear that, although electrolytes contain oppositely charged 
ions, these electrolytes may dissociate into their constituent ions 
which move about freely in the solvent medium. Arrhenius assumed 
an equilibrium between the free ions and the undissociated molecules 
of an electrolyte and sought to account for the observed properties 
of such a solution through displacement of this equilibrium. While his 
theory was successful in accounting for solutions of weak electrolytes, 
it led to contradictions in the case of strong electrolytes. These con- 
tradictions were not resolved until Debye introduced his theory of 
ion interaction some fifteen years ago. 


2. Ion interaction: thermodynamic properties. The properties of 
electrolytic solutions fall into two classes: (1) those of reversible and 
(2) those of irreversible processes. The first class includes all thermo- 
dynamic properties such as freezing point, boiling point, vapor pres- 
sure, and the various thermodynamic energy functions; the second 
includes all time-dependent phenomena, such as conductance, diffu- 
sion, and transference. 

According to Gibbs, equilibrium in a multiphase system requires 
that the thermodynamic potential of any component shall have the 
same value in the several phases. In systems in which reaction occurs 
and equilibrium exists among the reactants, the sum of the thermo- 
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dynamic potentials of the reactants must be equal to zero. It becomes 
necessary, therefore, to determine how the thermodynamic potential 
of an electrolyte depends upon the various thermodynamic variables. 

In a dilute solution of a substance between whose molecules no 
forces act (an ideal solute), the thermodynamic potential of the solute 
varies as the logarithm of the concentration; thus 


wi = RT log C +7, 


where yp; is the thermodynamic potential of the ideal solute, R the 
gas constant, J the absolute temperature, C the concentration in 
moles per liter, and J an integration constant. If forces act between 
the molecules of the solute, then the observed thermodynamic po- 
tential becomes 


uo = RT logC + RT log f +1, 


where RT log f is the energy due to the interacting forces between the 
molecules. It is evident that this energy is RT log f=po—wpi. 

If we have an electrolyte AB which is in equilibrium with its ions 
A+ and B~ according to the reaction equation A++B-=AB, then, 
according to the theory of Arrhenius, we should have the equilibrium 
equation 

Cat Ca 
Cas 


where K is the equilibrium constant. This equation is based on the 
assumption that no forces act between the ions, between the ions 
and the undissociated molecules AB, or between the molecules AB 
themselves. In the light of our preceding discussion, this assumption 
is not permissible, since the Coulomb forces, which vary as the inverse 
square of the distance, are constantly acting between the ions. The 
difficulty may be resolved by taking into account the energy due to 
the interaction between the ions. For a binary electrolyte, this en- 
ergy term is written in the form RT log f?, in which case our equi- 
librium equation assumes the simple form 


Cat X Cer 


2 = 
Cas 


If y is the fraction of the electrolyte dissociated into ions, then 
and Cas=C(1—vy). Hence, our equation takes the 
form 

Cy? 
(1) 
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If we know the value of f (the activity coefficient), it becomes possible 
to calculate the equilibrium constant K, provided y is known. 

Let us first consider the case of a completely dissociated electrolyte, 
(vy =1), at relatively low concentration; the positive and negative ions 
will be distributed throughout the solvent medium in a somewhat 
random manner due to their energy of thermal agitation. This dis- 
tribution of the ions, however, will not be the same as that of un- 
charged molecules between which no forces act. Due to the Coulomb 
forces acting between the ions, there is a greater probability that a 
positive charge will be found in the neighborhood of a negative charge 
than in the neighborhood of another positive charge, and the same is 
true of a negative charge in the neighborhood of a positive charge. 

Debye has succeeded in computing the charge distribution in the 
neighborhood of an ion. A detailed account of the method employed 
would be out of place here, but I may briefly present the physical 
picture and give the final result. Let us fix our thoughts for the mo- 
ment on a positive ion bearing a unit charge; it will, on the average 
(with respect to time) be surrounded by one unit of negative charge, 
distributed symmetrically with respect to the central ion, and the 
amount of charge in a shell of thickness dr, around the positive 
charge, will vary in a regular manner with the distance from the 
central charge. Actually, the charge in a shell of given thickness grad- 
ually increases from very low values, passes through a maximum, and 
then approaches zero at great distances. According to the Debye 
theory, the charge density falls off to 1/eth its value at a distance 
1/x=1.99X10-'!°(DT/C)/?, where D is the dielectric constant of the 
medium. As the concentration diminishes, the distance 1/x increases, 
and it follows that, in diluting a solution, work must be done against 
the forces acting between the central ion and its surrounding charge, 
or its ion atmosphere, as it is called. 

For a relatively dilute solution, the equation for the activity coeff- 
cient of the ions takes the simple form 


(2) — logio f = BC", 
where 
1.31 X 10° 
(DT)3!2 


In our equations above, we have assumed that each ion carries a 
unit charge. Similar relations are easily obtained for the case of ions 
with any number of charges. For the sake of simplicity, our subse- 
quent discussion will be restricted to uni-univalent salts, unless other- 
wise mentioned. 
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In Fig. 1 are plotted values of —log f against the square root of the 
(total ion) concentration for uni-univalent, uni-divalent, and uni-tri- 
valent salts. As may be seen, at lower concentrations the points ap- 
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Plot of —log f as a function of the square root of total ion concentration 
for salts in water. 


proximate closely to the straight lines required by theory. As the 
charge on one of the ions increases, the curves become steeper in ac- 
cord with theory. The equation has been tested with the very exten- 
sive experimental material relating to aqueous solutions and has been 
found to hold over a considerable range of concentration in the dilute 
region. At high concentrations, however, the deviations become large. 


3. The Debye-Hiickel-Onsager equation. Let us now turn to the 
irreversible phenomenon of the electrical conductance of a solution 
of an electrolyte. This is an important case, since the conductance is 
the one property of electrolytic solutions that may be measured with 
high precision over the widest range of conditions. It is known that 
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in the case of all strong electrolytes the equivalent conductance A 
(the specific conductance divided by the equivalent concentration) 
increases with decreasing concentration and approaches a limiting 
value. Arrhenius accounted for the gradually decreasing conductance 
at higher concentration by assuming an equilibrium between ions and 
undissociated molecules, the equilibrium being shifted with increas- 
ing concentration in the direction of a larger proportion of undisso- 
ciated molecules. Debye and Hiickel, however, showed that the 
decreased conductance at higher concentration is due to a decrease 
in the mobility of the ions themselves. The reason for this is not diffi- 
cult to understand. 

The force acting on an ion depends not only upon the external field 
under which it tends to move, but also upon the field due to the ion 
atmosphere. In the case of a stationary ion, the ion atmosphere is 
symmetrical with respect to the central charge, and the resultant 
force acting on the charge is zero. With a moving ion, however, the 
ion atmosphere becomes unsymmetrical with respect to the central 
charge and exerts upon this charge a force that counteracts the ap- 
plied field and thus tends to reduce the mobility of the ion. The asym- 
metry of the ion atmosphere of a moving charge is due to the fact 
that it requires a finite time to build up the ion atmosphere about a 
charge and, similarly, it requires a finite time for the ion atmosphere 
about a charge to become dissipated when the central charge disap- 
pears. This time is called the relaxation time and may be computed 
from known constants of the electrolyte and of the medium in which 
it is dissolved. Since the diameter of the ion atmosphere (1/k) di- 
minishes with increasing concentration, the asymmetry of the charge 
(in the ion atmosphere) exerts an increasingly great influence upon 
the moving ion as the concentration increases. 

Another factor needs to be taken into account, namely, what is 
known as the electrophoretic effect. Due to the charge of the ion at- 
mosphere, the water molecules in this atmosphere are carried in a 
direction opposite to that in which the ion moves. This results in a 
further diminution in the forward movement of the central charge. 

The influence of the ion atmosphere on the conductance of an elec- 
trolyte, as computed by Debye and Hiickel and Onsager, takes the 
comparatively simple form of the following equation: 


(3) A = y(Ao — a(;C)'/?), 


where a=8.18 X105/(DT)*/2A,+82/n(DT)*/?. Here 7 is the viscosity 
of the solution, and Ao is the limiting value which the equivalent con- 
ductance approaches as the concentration is diminished without 
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limit; in other words, it is the conductance of the ion when the in- 
fluence of the ion atmosphere becomes vanishingly small. 

In a solution of a strong electrolyte in a solvent of high dielectric 
constant (such as water), the electrolyte is completely dissociated, 
and y, in the above equation, becomes equal to unity. In the equa- 
tions originally developed by Debye and Hiickel, y was assumed to 
be unity, and the correctness of this assumption has been amply con- 
firmed by the experimental facts in the case of strong electrolytes. 
In the case of weak electrolytes, however, it is necessary to introduce 
the degree of dissociation y. 

How well the Debye-Hiickel-Onsager theory accounts for the con- 
ductance of solutions of electrolytes in water is illustrated in Fig. 2, 
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Plot of equivalent conductance as a function of square root of concentration 
for NaCl and KCI in water at 25°C. 
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where the equivalent conductance of aqueous solutions of potassium 
chloride and sodium chloride are plotted against the square root of 
concentration. As may be seen from the figure, the experimental 
points lié upon curves that approach definite limiting tangents, shown 
in the figure as broken lines. These tangents, as drawn, are deter- 
mined by the coefficients a and are, therefore, not arbitrary. It is 
evident that, while the Debye-Hiickel-Onsager equation gives a cor- 
rect interpretation of the facts in very dilute solutions, at higher 
concentrations, even as low as 0.001 N, the deviations become ap- 
preciable and become increasingly larger as the concentration in- 
creases. 


4. Ion association. As stated above, for solutions of weak acids and 
bases it is necessary to assume incomplete dissociation, even in aque- 
ous solution. In other words, for these electrolytes the ions interact 
at short distances to form neutral molecules. There is reason for be- 
lieving that in the case of acids and bases something more is involved 
than merely the Coulomb forces acting between the ions. For normal 
electrolytes (which are strong electrolytes in water) dissolved in sol- 
vents of lower dielectric constant, we have good reason to believe that 
association of ions occurs due to Coulombic interaction at short dis- 
tances. Bjerrum considered this problem from the theoretical point 
of view as early as 1926. He showed that in solvents of lower dielec- 
tric constant the distribution of ions of one charge in the neighbor- 
hood of an ion of opposite charge takes the form shown in Fig. 3. 
(The figure is due to Dr. Fuoss.) Here there is a high concentration 
of ions within one or two ion-diameters of the central ion, then a 
minimum in the distribution curve, and thereafter a maximum simi- 
lar to that in the ion atmosphere of strong electrolytes in water. If it 
is assumed that an ion lying near to the central ion takes no part in 
the conduction process, while an ion at a distance exerts its normal 
influence in the ion atmosphere, and that there is an equilibrium be- 
tween the two types according to the law of mass action, then the 
dissociation constant of the electrolyte may be related to the dielec- 
tric constant of the solvent medium and the sum of the radii of the 
two ions. The complete expression is somewhat complex, but a first 
approximation takes the form 


1 4xnN aDkT e? 
a’ exp ( 
K 1000 e aDkT 


where N is Avogadro’s number, & is Boltzmann’s constant, e is the 
unit charge, and a is the distance between the centers of charge of the 
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two ions when in contact. It is seen from this equation that the dis- 
sociation constant K is closely related to the dielectric constant D of 
the medium, on the one hand, and to a, the distance of closest ap- 
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Distribution of charge around central ion in solvent of dielectric constant 
D=20. a=5.57 X10-®. C=0.5, 1, 2, 310-3. 


proach of the ions, on the other. It is evident that the value of K 
diminishes rapidly with decreasing value of the distance a, or with 
decreasing value of the dielectric constant D. 

In Fig. 4 is shown the curve connecting —log K with the dielectric 
constant. The curve is based on the value a=6.4 X10-* cm. and the 
circles represent experimental points for solutions of tetratsoamylam- 
monium nitrate in mixtures of dioxane and water. As may be seen 
from the figure, the dissociation constant reaches values approaching 
unity for a value of D in the neighborhood of forty. For this and 
higher values of D, the electrolyte is completely dissociated. For 
D=2.4 the dissociation constant has fallen to 2107". 
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We should expect that all electrolytes in solvents of lower dielectric 
constant would exhibit marked association of their ions, and that the 
association would be the greater the smaller the ions. This has been 
amply confirmed by the work of Dr. Fuoss and others in the Brown 
laboratories. 
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Plot of —log K as function of log D. a=6.4X10-*. 


In examining experimental conductance data in the light of the 
above theoretical considerations, it is necessary to solve simultane- 
ously equations (1), (2), and (3) above. While an explicit solution is 
not possible, Dr. Fuoss has developed ready methods of approxima- 
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tion that render it possible to evaluate the experimental results in 
terms of the theory. Essentially, what is necessary, according to the 
present theory, is to take into account the equilibrium between free 
and associated ions according to the law of mass action, and the long 
range interaction of the ions between one another in terms of the ion- 
atmosphere effect. How well the theory accounts for the experimental 
values is shown in Fig. 5, where are plotted values of AoF/A against 
those of CAf?/F, where f is the activity coefficient of the ions and F 
is a function readily computed from conductance data (Fuoss, Jour- 
nal of the American Chemical Society, vol. 57 (1935), p. 488). 
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As may be seen from the figure, the experimental points closely ap- 
proximate the linear relation required according to theory. At higher 
concentrations, however, the deviation becomes increasingly larger. 
Again, theory supplies a satisfactory limiting solution. 

For solutions of electrolytes having spherical ions, the dissociation 
constant is determined by a, the sum of the ion radii, and D, the di- 
electric constant of the solvent medium. If K is known, a may be 
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computed therefrom. The values of a so obtained are in reasonable 
accord with atomic dimensions as determined by other means. With 
electrolytes having complex or unsymmetrical ions, the distance a de- 
pends upon the relative position of the ions and the simple theory 
(based on spherical ions) becomes an approximation. Nevertheless, 
the experimental values are fairly well accounted for even in this case. 

In Table 1 are given values of Ao, K, and a for a number of elec- 
trolytes in liquid ammonia. 


TABLE 1 


Conductance of salts in liquid ammonia 


Salt Ao K aX 108 
LiNO; 277 36.5 5.04 
NaNO; 295.1 28.8 4.46 
KNO; 329.0 14.7 3.57 
RbNO; 341.8 11.4 3.24 
CsNO; 333.5 9.66 3.21 


As appears from the table, electrolytes with larger ions have lower 
conductance values and larger dissociation constants. The ion diame- 
ters are reasonably in accord with atomic dimension. The reversal of 
size in the case of the alkali metal ions is doubtless due to solvation. 


5. Ion-dipole interactions. Before proceeding further, it will be 
helpful if we examine the conductance curves of an electrolyte as a 
function of concentration in solvent media over a wide range of the 
dielectric constant. Conductance curves are shown in Fig. 6 for solu- 
tions of tetraisoamylammonium nitrate in mixtures of dioxane and 
water whose dialectric constants vary from 2.2 (dioxane) to 78 
(water). Logarithms of concentration and of equivalent conductance 
are plotted in order that all values may be represented in a single 
figure. The concentration varies from approximately 0.1N to 
110-5N, and the conductance from 100 for water to 310-5 
for dioxane. The first two curves (D=78 and 37) are normal, ap- 
proximating the Debye-Hiickel-Onsager equation. The third curve 
(D=12) shows a large conductance decrease; here ion association is 
beginning to occur. With the fourth curve (D=9), a minimum ap- 
pears in the curve at about 0.03 N. As the dielectric constant dimin- 
ishes, the conductance values fall and the minimum shifts to lower 
concentrations, reaching 2.5 10-5 in pure dioxane. For dilute solu- 
tions, in mixtures having values of D less than 6, the curves approach 
limiting tangents with a common slope of —1/2. This corresponds to 
a mass action equilibrium between free and associated ions of a 
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largely associated electrolyte (y small). The complexity of the lower- 
most curves at concentrations above that of the minimum point are 
of much interest. I shall refer to them again. 
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Plot of log A as function of log C for solutions of tetratsoamylammonium 
nitrate in water-dioxane mixtures of varying dielectric constants. 


The minima in the conductance curves are due to interaction be- 
tween the simple ions and the associated ions or ion pairs. Thus, if 
we have an electrolyte AB, dissociating into the ions A+ and B-, the 
ions A+ and B- will interact with the ion pairs AB. An equilibrium 
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exists between them according to the equations AB+At=ABAt 
and AB+B-= BAB-. This interaction has been treated by Dr. Fuoss 
in much the same manner as that between ions and ions. It leads to 
the following simple equation for the conductance as a function of 
concentration: 


X Ki/2 A 
(4) AC¥2g(C) = + — (1 


Here g(C) is a correction factor taking into account the activity and 
the mobility terms of the ions; Xo is the limiting conductance of the 
triple ions ABA+ and BAB-, assumed to be equal; and is the dis- 
sociation constant of the triple ion equilibria. Since A/Ao is small in 
comparison with unity, a plot of the AC'/*g(C) against C should yield 
a linear relation. A plot of this kind is shown in Fig. 7 for solutions 
of tetrabutylammonium picrate in anisole. As may be seen from the 
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Plot of conductance function for tetrabutylammonium picrate in 
anisole showing effect of triple-ion equilibrium. 


figure, the experimental points lie upon the required straight line well 
within the limit of experimental error. At higher concentrations, 
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slightly above the concentration at the minimum point, deviations 
become appreciable and increase with increasing concentration. 


6. The effect of potential. An ion moving under the action of a 
constant applied field is retarded owing to the asymmetry of its ion 
atmosphere. On the one hand, as we have seen, a potential develops 
opposing the potential of the field, and, on the other hand, the elec- 
trophoretic effect of the ion atmosphere tends to move the solvent in 
a direction opposite to that of the motion of the ion. Since the time 
required to build up the ion atmosphere is finite, the condition arises, 
when the applied field becomes sufficiently strong, that the ion moves 
away from its atmosphere more rapidly than it (the ion atmosphere) 
can be built up. In the limiting case, the ion moves quite free of an ion 
atmosphere and it is in the same condition as it would be if the solu- 
tion were infinitely dilute. Since the relaxa ion time of the ion atmos- 
phere varies inversely as the concentration, this field effect is greater 
(for a given potential) the lower the concentration of the solution. 
This effect was discovered by Wien and has been studied extensively 
by him and his collaborators; the results are in accord with theory. 
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Illustrating the Wien potential effect for solutions of strong electrolytes in water. 


The effect of potential on the conductance of aqueous solutions of 
strong electrolytes is illustrated in Fig. 8, where applied potentials 
are plotted as abscissas against the increase of conductance due to the 
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field divided by the maximum increase possible due to an infinite 
field. Since the relaxation time is an inverse function of concentration, 
the conductance change on applying the field reaches its limiting 
value much more rapidly at low concentrations than at high concen- 
trations. 

In the case of weak electrolytes, another effect (due to the field) 
appears. Owing to the field which tends to separate the oppositely 
charged ions in their ion pairs, the dissociation of the electrolyte is 
increased as the field increases. This effect has also been studied by 
Wien, and the theory has been developed by Onsager. 

In Fig. 9 are shown curves for the percentage conductance change 
of several weak acids in aqueous solution. The experimental points 
are indicated and the straight lines, as drawn, represent the course 
of the curves required by theory. 
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Illustrating the Wien potential effect for solutions of weak electrolytes in water. 


7. Frequency effect. The effect of potential on conductance led 
Debye and Falkenhagen to suspect an effect due to the frequency of 
the applied field. They developed a theory which was later confirmed 
experimentally. The effect is due to the fact that a finite time is re- 
quired for the asymmetry of the ion atmosphere around a moving 
ion to disappear when the field is removed. With a slowly alternating 
field the effect is inappreciable, but when the period of the field ap- 
proaches the relaxation time of the ion atmosphere, the asymmetry 
effect of the ion atmosphere on conductance diminishes and disap- 
pears at sufficiently high frequencies because the displacement of the 
central ion in its atmosphere becomes inappreciable in the period of 
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an oscillation. The electrophoretic effect due to the ion atmosphere 
itself remains. 

This conductance change over that of the conductance in a sta- 
tionary field decreases with decreasing concentration. It likewise de- 
creases with decreasing number of charges on the central ion. If the 
conductance of an electrolyte is measured at a given frequency com- 
parable with the relaxation time of the ions at intermediate concen- 
trations, then, at very low concentrations where the radius of the ion 
atmosphere is very large, the frequency effect will be zero, since the 
ions will already have a conductance equal to that of their limiting 
value for an infinitely dilute solution. With increasing concentration, 
the conductance change due to frequency increases, passes through a 
maximum, and thereafter decreases. The reason for the decrease is 
that the frequency of the applied field is now below that required in 
order to prevent the building up of the asymmetry of the ion atmos- 
phere. 

In Fig. 10 are shown curves for the percentage change in the con- 
ductance of magnesium sulfate and manganese sulfate in aqueous 
solution for a field frequency corresponding to a wave length of one 
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Illustrating the frequency effect for solutions of salts in water. 


meter. For magnesium sulfate, the curve approximates closely to the 
theoretical value over the entire range of concentration; for man- 
ganese sulfate, a large deviation occurs at the higher concentrations. 
In the case of simpler salts, such as uni-univalent salts, the corre- 
spondence between theory and experiment is even better than in that 
of the two salts above. 
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8. Ion dipoles. In measuring the conductance of solutions, we ob- 
tain information relating directly to the ions, but our knowledge with 
respect to the associated ions or ion pairs depends upon inference. 
It is possible, however, to carry out measurements in which the prop- 
erty measured depends directly upon the ion pairs and not upon the 
ions. 

The dielectric constant of a vacuum is unity; that of all material 
media is greater than unity. This is due to the fact that the charges 
in the atoms of molecules of the medium suffer displacement under 
the action of an impressed potential. With non-polar molecules, this 
displacement of charge is extremely small, but if the centers of charge 
of positive and negative electricity in a molecule are not coincident 
(polar molecules), then the molecule has an electric moment, and it 
will tend to rotate under the action of the field so as to align the axis 
of its electrical dipole with the axis of the field. The presence of such 
dipole molecules in a solvent medium tends to increase the dielectric 
constant of the solution over that of the pure medium. 

If the charges in the dipole are actually separated, as we should 
expect them to be in the case of associated ions, or ion dipoles, as we 
may call them, then we should expect that the effect of these dipoles 
upon the dielectric constant would be very large. The extent of the 
displacement of charge is stated in terms of the molecular polariza- 
tion of the dissolved dipoles, and from the polarization we may com- 
pute the electric moment (charge times distance) of the dipoles. If, 
in the dipole molecules, the charges are actually separated, as they 
are in the case of ionic dipoles, then by dividing the electric moment 
by the unit charge on the ions, we obtain therefrom the distance of 
separation of the charges. 

In Fig. 11 are shown the polarization curves of a number of elec- 
trolytes as derived from dielectric constant measurements of their 
solutions in benzene. These polarization curves are in part linear and 
slightly inclined, and in part they show a much higher degree of com- 
plexity. The inclination of the linear curves is probably due to the 
association of the dipoles to form quadrupoles. The rapid diminution 
of the polarization in the case of other electrolytes is doubtless 
due to a high degree of association of the dipoles to more complex 
molecules. 

What is striking, however, is the fact that the values of the polari- 
zation for electrolytes dissolved in benzene are approximately ten 
times as great as those of ordinary polar molecules in the same sol- 
vent. These high values of the polarization indicate that the charges 
are separated by considerable distances. Values of the polarization, of 
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Fic. 11 
Polarization curves of electrolytes in benzene solution. 


the dipole moments, and of the distance between charges, as com- 
puted from the dipole moments, are given in Table 2. 


TABLE 2 
Electric moments of ion dipoles in benzene 
Salt Polarization Moment X 1018 aX108 
- Ac 2690 11.2 2.34 
Bu,N- Br 3250 12.3 2.58 
Bu;NH - Pi 3770 13.2 2.78 
Am;NH- Pi 3830 2.80 
Am,N-SCN 5050 15.4 
Bu,N - Pi 6740 17.8 3.74 
Am,N- Pi 7090 18.3 3.83 
Am,N-HOBPh; 8200 19.6 4.12 


Bu=butyl, Am =amyl, Ph=phenyl, Ac=acetate, Pi=picrate. 
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On inspection of the table, it will be seen that the distance between 
charges varies from 2.3 to 4.1 Angstrém units (A=1X10-* cm.). 
These values are in reasonably good accord with what we know of 
molecular and atomic dimensions. It is particularly interesting to 
note that as the ions become more complex, and therefore larger, the 
distance between the ions increases. 


9. Molecular weights. The decrease in the molecular polarization 
of electrolytes in benzene with increasing concentration indicates the 
association of dipoles to quadrupoles and perhaps more complex 
structures. Dielectric constant measurements, however, can give us 
no direct evidence as to the existence of these structures. Some fur- 
ther light is thrown on the phenomenon of association by measure- 
ments of the molecular weights of electrolytes in benzene. While the 
force between ions and ions falls off as the inverse square of the dis- 
tance, and that between ions and dipoles as the inverse cube of the 
distance, the force between dipoles and dipoles falls off as the inverse 
fourth power of tre distance. We should not, therefore, expect to find 
a high degree of dipole association in solvents of higher dielectric con- 
stant; but in solvents of very low dielectric constant, such as benzene, 
we should expect association of dipoles to occur. Indeed, we might 
look for association complexes of higher order, including charged as 
well as uncharged complexes. If dipole association occurs, measure- 
ment of the molecular weight of the dissolved electrolyte should indi- 
cate, approximately, its extent. In this connection, it may be noted 
that the concentration of ions in solutions of electrolytes in benzene 
is so low that the ion atmosphere effect may be neglected. 

In Fig. 12 are shown curves in which are plotted values of the 
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Freezing point curves of solutions of electrolytes in benzene. 
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ratio At/At) of the observed freezing point lowering to the lowering 
of an ideal solute having the same (stoichiometric) concentration. 

It will be observed that in the case of all electrolytes the curves lie 
below that of an ideal solute, indicating a measurable amount of asso- 
ciation. An analysis of the curve for tributyl-ammonium picrate (the 
least associated) indicates that there exists an equilibrium between 
dipoles and quadrupoles in accordance with the law of mass action. 
In the case of other electrolytes, association takes place to a larger 
extent, and the observed values cannot be accounted for on the as- 
sumption of a simple equilibrium between dipoles and quadrupoles. 
Evidently, in these cases, higher complexes than quadrupoles are 
formed. This is particularly true of tetratsoamylammonium thiocy- 
anate, where the mean molecular weight of the solute at a concentra- 
tion of 0.017 N is eleven times its formula weight. There can be no 
doubt but that very complex aggregates are formed in solutions of 
this electrolyte. 

The results of these molecular weight determinations are in accord 
with those of dielectric constant measurements. It will be recalled 
that the molecular polarization of tetratsoamylammonium thiocya- 
nate falls off very sharply with increasing concentration. The results 
are also in accord with the results of conductance measurements. In 
Fig. 6, for example, for tetratsoamylammonium nitrate, the conduc- 
tance curve for low values of D shows several inflection points at 
concentrations above the minimum; there is a certain region where 
the conductance increases very slowly with concentration. A similar 
curve in the case of tetratsoamylammonium thiocyanate shows an 
even more marked effect; there is a considerable region where there 
is scarcely any increase in conductance as the concentration is in- 
creased. Presumably, what is happening here is the formation of a 
large proportion of highly complex neutral structures which, while 
taking no part in conductance, tend to reduce the normal increase 
of ions with increasing concentration. 

Since in all cases the conductance ultimately rises very rapidly 
(at high concentrations), it follows that in these regions there must 
be a rapid increase in the number of ions. Seemingly, there is an 
equilibrium among a great variety of complex ions and neutral mole- 
cules which, with increasing concentration, is shifted in the direction 
of an increased number of charged structures, that is, structures in 
which there is a preponderance of either positive or negative ions. 


10. Conclusion. The preceding discussion will have shown that the 
theory based on the assumption of Coulombic interaction between 
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ions gives an uncommonly satisfactory account of the properties of 
electrolytic solutions in the limiting case of dilute solutions. It must 
not be forgotten, however, that the theory is strictly a limiting the- 
ory, and that as yet we have no adequate theory of more concentrated 
solutions. While a large mass of data is available with respect to con- 
centrated solutions, an adequate theory is lacking. That the present 
theory may be expanded in some simple fashion so as to account for 
the more concentrated solutions seems doubtful. It appears more 
probable that an entirely new method of attack will have to be de- 
vised. 

Actually, only a beginning has been made in the solution of the 
problem of electrolytes. The nature of the ions themselves remains 
unknown. We know that ions in water are hydrated, but the extent 
of the hydration is not known. There is evidence indicating that ions 
in media other than water are solvated, but again the extent of 
solvation remains unknown. 

It is known that the conductance of an electrolyte in a given 
medium is roughly proportional to the fluidity of the medium, but 
the precise relation between conductance and viscosity remains un- 
determined. So also the relation between the conductance phe- 
nomenon and temperature remains very obscure. Indeed, with the 
exception of the measurements carried out by A. A. Noyes and his 
co-workers thirty years ago, we have practically no data with re- 
spect to aqueous solutions at temperatures above 100°, and even at 
that temperature the data are meagre and uncertain. In solvents 
other than water, there is almost a complete lack of data over any 
considerable temperature range. 

The generality of electrolytic phenomena is still not recognized by 
many investigators. Practically all liquid media are capable of dis- 
solving suitable electrolytes and yielding solutions that conduct the 
electric current. Much work remains to be done with solvents of very 
low dielectric constant. Data with respect to solutions in solvents of 
dielectric constant higher than that of water are also lacking. 

Many substances which are not themselves electrolytes in the pure 
state yield electrolytic solutions when dissolved in suitable solvents. 
Such is the case with the acids and many bases in water. Even today 
we have no sound theory to account for the influence of constitution 
upon the strength of acids and bases. The properties of acids and 
bases in solvents other than water remain uninvestigated. So also 
the properties of the weaker salts in solvents other than water remain 
unknown. The investigation of the remarkable solutions of the alkali 
metals in liquid ammonia has been carried only far enough to show 
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that these solutions form a link between metallic and electrolytic con- 
ductors. These solutions offer a tempting field for further research. 

Electrolytic solutions find applications in many diverse fields of 
science, particularly in the biological sciences. Biological systems, in 
general, are electrolytic systems and many biological phenomena di- 
rectly involve electrolytes. It does not seem probable that biological 
processes may be accounted for merely on the basis of what is known 
of aqueous solutions. Cell structures consist of anything but water 
and yet these cell structures are of primary importance in connection 
with many electrolytic phenomena that accompany biological change. 
The transmission of nerve impulses is recognized to be an electrical 
phenomenon. Clearly, the electric impulse in this case is nota metallic 
one; the process must, necessarily, be an electrolytic one. An under- 
standing of the mechanism of such phenomena will be had only when 
the nature of electrolytic processes in cells is more thoroughly under- 
stood. 
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TORSION OF REGIONS BOUNDED BY 
CIRCULAR ARCS* 


I. S. SOKOLNIKOFF AND E. S. SOKOLNIKOFF 


1. Introduction. N. MuscheliSvili discussed recentlyt some elegant 
general methods of solution of a class of problems occurring in the 
theory of elasticity. These methods will bear closer scrutiny as the 
present authors have convinced themselves by obtaining, with re- 
markable ease, solutions of St. Venant’s torsion problem for some 
regions bounded by circular arcs. A number of problems in harmonic 
and biharmonic analysis, recently discussed by other investigators, 
can be treated with greater simplicity by these methods. 

Muschelisvili makes use of the fact that a function of a complex 
variable {, analytic for |{| <1, whose real part on |{| =1 assumes 
continuous real values u=x(@), (0<0<2z), is given by 

1 u(6)do 
(1) = — 


+ const., 


where =e, and ¢ is any point interior to the unit circle y. It is 
readily shown that (1) is entirely equivalent to the formula of 
Schwarz. 

Consider a continuous, simple closed curve C in the complex z 
plane and let z=w({) map the region bounded by C conformally on 
the unit circle y in the ¢ plane. Let F(z) be analytic in the interior of 
C, and let the value of the real part of F(z) on the boundary of C 
be u(x, y). 

Substituting 
w(t) + _ w(t) — 


2 7 2i 


in u(x, y) gives 
1\ 
u(x, y) a) —); 


where bars denote the conjugate values. Hence, from (1), 


* Presented to the Society, December 30, 1937. 

¢ N. Muscheli8vili, Rendiconti delle Reale Accademia dei Lincei, vol. 9 (1929), 
pp. 295-300; Mathematische Annalen, vol. 107 (1932), pp. 282-312; Zeitschrift fiir 
Angewandte Mathematik und Mechanik, vol. 13 (1933), pp. 264-282. 
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1 

1 («. -) do 
(2) F(z) = F[w(¢)] = —f + const., 

Tid, 
and the problem of Dirichlet for the region bounded by C is solved. 

It may be remarked that the procedure just outlined differs from 

that formerly used by one of the present authors* in that the integral 
of Poisson employed by him is replaced by a modified integral of 
Schwarz. The latter frequently admits of easy evaluation, in closed 
form, as is illustrated below in the treatment of a torsion problem. 


2. Torsion problem. Consider a region R of the complex z plane 
bounded by two circular arcs C; and C2 making an angle a0 at 
their points of intersection z =z, and z =z. The function effecting con- 
formal mapping of the region R on the unit circle in the complex 
¢ plane is 


Ci(z — 23)7/* + (2 — 


(3) 


The choice of the constants C, z:, and z is uniquely determined by 
the geometrical configuration and the scale used in mapping the re- 
gion R on the unit circle. For a=2/n, where 7 is an integer, the map- 
ping function becomes rational, which fact greatly simplifies the eval- 
uation of (2). To illustrate the procedure it will suffice to consider 
two important special cases when the region R is one of the following: 
(a) a lune formed by two circular arcs of equal radius and inter- 
secting at right angles; 
(b) a semi-circle. 
The first of these regions suggests the profile of an airplane wing. 
It is easily verified that the mapping functions appropriate to the re- 
gions defined in (a) and (b) are, respectively, 


22 


(4) 


and 
+ 1)? — i(z — 
(z + 1)? + i(z — 1)? 


In obtaining (4) and (5) the constants 2, and z, were chosen as —1 
and +1, which can be done without loss of generality inasmuch as 


(5) 


* I. S. Sokolnikoff, Transactions of this Society, vol. 33 (1931), pp. 719-732. 


386 I. S. SOKOLNIKOFF AND E. S. SOKOLNIKOFF [June 


this choice merely fixes the scale for the radii of arcs bounding the 
regions. The constant C was chosen so as to make the point z=0 in 
(4) correspond to ¢ =0, and the point z=0 in (5) correspond to ¢ = —1. 

The transformations (4) and (5) permit of easy solution, in closed 
forms, of the torsion problem for prisms whose cross sections are the 
regions defined in (a) and (b). The solution of the second of these 
problems is included in the papers of Greenhill,* who attacked the 
problem in an entirely different way. The solution of the first problem 
appears to be lacking, and it will be given next. 

Since the conjugate V(x, y) of the torsion function ®(x, y) assumes 
on the boundary of the cross section of the prism the value 


W = 3(x? + y*) + const. 
= 423 + const. = }w(c)w(c) + const., 
it follows from the foregoing that the real part of the analytic function 
1 w(c)w(o 
(6) =— + const. 
o 


gives the desired function V(x, y). From (4) it follows that w(c) 
= [1—(1—o*)?]/o, where the appropriate branch of the square root is 
determined from the observation that the imaginary part of z is posi- 
tive whenever @ in ¢ =e* lies between 0 and 7. 

Noting that on |¢| =1 we have ¢=1/e, we deduce the relation 


w(o)w(e) = [1 — (1 — E 


Substituting this expression in the integrand of (6) and evaluating 
the resulting elementary integrals, we obtainf 


i+? 
or if we return to the z-plane, 
s*— 1 (1 — s*)? i-s 


F(z) = log + const. 


2+1 +2) 1+2 


* A. G. Greenhill, Messenger of Mathematics, vol. 8 (1879), p. 89; vol. 10 (1881), 
p. 83. 

+ The details of integration are of little interest and hence are omitted. In evaluat- 
ing the integrals involving square roots, it is advantageous to make a substitution 
o=(1+7f)/(1—7), so that the resulting integrals are evaluated over the real axis 
of the ¢ plane. 
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The real part of F(z), which is the desired solution, is 
9}, 
where 
[(1 — x? — y?)?4 4y?]1/2 — 2y 
= and ¢ = - 
(1 + x)? + 1— x? — y’? 
A simple calculation shows that, on the boundary of the lune 
formed by x?+(y+1)?= 2, V(x, y) reduces to 


+ y? 
v= pine + const., 


as it should. 
An analogous calculation* gives for the region bounded by a semi- 
circle 


ri(i + 9) 


f) = 


2 (1 i(1 — 

log 

mi (t+ $)? 1+¢ 
or, if we return to the z-plane with the aid of (5), 


2(27 + 1 2? — 1)? 
( ) log + const., 
z 2? 1+ 2) 


+ const., 


1 
F(z) = + 
2rt 


which agrees with the result of Greenhill in the form given by Love.t 

A similar procedure can be employed in treating St. Venant’s flex- 
ure problem, the solution of which, for the regions in question, does 
not appear in the literature.t{ 


THE UNIVERSITY OF WISCONSIN 


* The integrals appearing in this case are precisely of the same type as those oc- 
curring in the preceding problem. They are evaluated with the aid of the elementary 
theorems on residues, and by means of the substitution given in the preceding foot- 
note. 

+ A. E. H. Love, A Treatise on the Mathematical Theory of Elasticity, 4th ed., 1927, 
p. 319. 

t See, however, M. Seegar and K. Pearson, Proceedings of the Royal Society, 
Section A, vol. 96 (1920), pp. 211-232. 
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A MULTIPLIER RULE IN ABSTRACT SPACES* 
HERMAN H. GOLDSTINE 


The Lagrange multiplier rule has been generalized to certain non- 
calculus of variations problems by Graves,f Hahn,f{ and the author.§ 
Moreover a very general problem was formulated by L. A. Luster- 
nik.|| However his work seems to rest upon a theorem which is stated 
without proof and which the author is unable to verify. There are 
also certain other difficulties with his proof. The problem herein pre- 
sented is so formulated that the problem of Bolza in the calculus of 
variations, the problems treated by Graves, Hahn, and the author, 
and the problem of minimizing a functional defined upon an arbitrary 
Banach space subject to very general numerically-valued side condi- 
tions, and numerous other examples are included as special cases. 

The proof proceeds along lines which are essentially generalizations 
of the methods of the calculus of variations. This demonstration is 
made possible by the very powerful implicit function theorems of 
Hildebrandt and Graves] and yields analogs of the transversality 
condition and of the Euler-Lagrange equations. Some instances which 
explain the number of linear spaces involved are given in the conclud- 
ing section. 


1. Definitions and assumptions. In order to obtain our statement 
and proof of the multiplier rule it is convenient to refer to five normed 
linear spaces M, N, ¥, U, BV, of which M, U, and B are assumed to 
be complete. It is further supposed that Mo, Mo, and Xo are regions 
of their respective spaces. We shall then be concerned with the fol- 
lowing basis: 


* Presented to the Society, September 10, 1937. 

+ A transformation of the problem of Lagrange in the calculus of variations, Trans- 
actions of this Society, vol. 35 (1933), pp. 675-682. 

t Ueber die Lagrange’ sche Multiplikatorenmethode, Sitzungsberichte der Akademie 
Wien, vol. 131 (1922), pp. 531-550. 

§ The minima of functionals with associated side conditions, Duke Mathematical 
Journal, vol. 3 (1937), pp. 418-425. 

|| Sur les extrémes relatifs des fonctionnelles (in Russian), Recueil Mathématique 
de la Société Mathématique de Moscou, vol. 41 (1934), pp. 390-401. 

{ See L. M. Graves, Topics in the functional calculus, this Bulletin, vol. 31 (1935), 
pp. 641-662; Riemann integration and Taylor's theorem in general analysis; Implicit 
functions and differential equations in general analysis; and T. H. Hildebrandt and 
L. M. Graves, Implicit functions and their differentials in general analysis. The last 
three articles cited appear in the Transactions of this Society, vol. 29 (1927), pp. 163- 
177, 514-552, and 127-153, respectively. In the sequel these papers will be referred 
to as papers I to IV, respectively. 
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A. f and yx, (k=1, - --, g), where g is an arbitrary integer, are 
real-valued functions defined on the composite* (Xo, Dto) of ¥o and 
Mo, each of which has a first differential at some fixed point (Zo, fo). 
Moreover the functionals ¥;, (k=1,---, g), vanish at this point 
fo). 

B. u=¢i(xo, Mo, Yo) is an operation of class C’ defined on 
(Xo, Mo, No), and having its functional values in U1; and vp=n(xo, po) 
is defined over (Xo, Mo), has its values in Jo, is of class C’, and is such 
that Zo, (Zo, Zo) the zero element in the space U. (An- 
alogous notations will be used for the zero points in the various spaces 
treated.) 

C. There is an operation v= @e(Xo, uo, vo) of class C’, which is de- 
fined on a neighborhood of the point fo, where = (Zo, fo), 
and has its functional values in %, such that the linear and continu- 
ous function of u, (dudi(xo, Yo; Mo, Vo; has a linear 
continuous reciprocal at (xo, Mo, ¥o) = Ho, Fo). 

D. The functional f has a minimum at (%o, Zo) in the class of points 
in (Xo, Mo) which satisfy the equations 


$1[ xo, Ho, 2(Xo, po) | = 0,, 
wo) = 0, 


In order to simplify the notation the element i of B will be used 
to designate $2(%o, fo, Fo); Pa(x, u, v) to indicate d.(%o, fio, Fo; X, 
(a=1, 2); and N(x, to connote dn(Zo, fio; x, 


2. The imbedding theorem. It is of paramount importance for our 
purposes to show that the minimizing element (%o, fio) is a limit point 
of elements (xo, 40) which satisfy the equation 


Ho, | = Ou. 


To show that the element (%o, fo) is so imbedded it is convenient 
first to demonstrate several lemmas. 
LemMA 2.1. In a sufficiently small neighborhood of jio the equations 
(2.1) $1(%0, Ho, = 
$2(%0, Ho, Yo) = 
have, for each (xo, vo, u, v) in a neighborhood of (Zo, Fo, Ou, 3), @ unique 
solution, to = M,(Xxo, vo, u, v), of class C’. 


This result follows at once from the implicit function theorems 


* See paper IV for the definitions of terms such as composite, neighborhood, 
region, differential. 
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mentioned above,* because of assumptions C and the fact that 
(uo, Xo, Yo, ¥) = (flo, Zo, Fo, Ou, is an initial solution of these equa- 
tions. 


COROLLARY 2.2. For each (x, v, U, V) the equations 
#,(x, u,v) = U, 
$.(x, u,v) = V 
have a unique continuous solution p= M2(x, v, U, V) which is equal to 
dM,(Zo, Fo, Ou, 0; x, v, U, 


The existence and unicity of the solution is an immediate conse- 
quence of assumption C; and the last assertion in the corollary can be 
established by setting x9 = Zo+ax, vo=Fo+av, u=aU, and 
bo = M,(xo, vo, u, v) in equations (2.1) and differentiating the resulting 
identities in a at a=0. 

Then to obtain the desired imbedding theorems we make a final 
assumption. 

E. Either the function M,, mentioned above, is independent of vo, 
or the transformation 

p—d,M, [Zo, Bo, Ou, 5; dyn(Xo, Ho; 


has a linear continuous reciprocal. 

It follows readily from Lemma 2.1 that if M, is independent of 
vo and if M(xo, u, v) is defined to be Mi (xo, vo, u, v), then the functions 
M(xo, u,v), vol, b2[x0, M(xo, u, v), vo] are independent of vo, 
at least for (xo, vo, u, v) sufficiently near to (%o, Vo, Ou, 0). 


THEOREM 1. In a sufficiently small neighborhood of fio the equations 
Ko, n(Xo, po) | 
$2[-o, Ho, Ho) | 


possess a unique solution po= M(xo, u, v) of class C’ for (xo, u, v) in 
some neighborhood of Ou, 3). 


From what has been said above it suffices to consider the case 
where M, is not independent of vp. Hence the equation 


Mo = M,[xo, n(Xo, Mo), u, v| 


has a solution po = M(xo, u, v) which is effective in the theorem. 
Similarly a consideration of the equation 


t= dM, Ou, 0; N(x, U, 


* See paper I, p. 65, and paper IV, p. 150. 
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yields the following corollary: 
CoROLLARY 2.3. For each (x, U, V) the equations 
u, N(x, u)] = U, 
u, N(x, u)] = V 


have a unique continuous solution p= M;(x, U, V)=dM(%p, 3; 


COROLLARY 2.4. Let s be an arbitrary integer, (x;,---,x,) @ set 
of elements of X, and (u1,---, ws) @ set of elements of M satisfying 
[x, N(x, u)|=0,. Then there are functions xo(b) and o(b) of class 
C’ for values of b=(bi, - - - , b,) sufficiently near to (0, - - - , 0), which 
satisfy the equation 


$1[%0, Ho, = Ou, 
x0(0) =o, = fo, and = ms for 1 SiSs. 


Define V; to be ®,[x;, Bi, N(xi, |, (@=1,---, s), xo(b) 
b,) = +D;xi, and Ho(d) to be M [xo(b), I, 
where a repeated index indicates a sum. Then the first part of the 
corollary is obvious, and the latter part is an easy consequence of 
corollary 2.3. 


3. The multiplier rule. The usual mode of proof utilized in the 
calculus of variations suffices to establish the following lemma: 


3.1. There are constants lo, G1, Cg, not all zero, such 
that 


g(x, = lod f (Zo, Ho; x, u) + Zo; x. uw) = O 
for every (x, u) satisfying the equation 
,([x, B, N(x, 0.. 
To simplify the succeeding demonstration we shall define g(x) 
and ge(u) to be lod.f(%o, Ho; x) (Xo, Ho; x) and f(Zo, Ho; m) 
Zo; u), respectively. It is then obvious that g(x, = g:(x) 


+o(u) and that g; and g are linear continuous functionals on ¥ and 
M, respectively. 


LEMMA 3.2. The functions ®,[x, M;(x, U, V)], N[x, M;(x, U, V)|> 
and g|x, M(x, U, V)]| are independent of x and V; the latter is a linear 
continuous functional J on U1; and for every x 


(3.1) gi(x) + go[M3(x, 0., 0,)] = 0. 
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It follows at once from Corollary 2.3 that ®, is independent of x 
and V. Hence we have ®,[#, z, N(z, where #=x1—x2, and 
a= M;(2, 0., Vi— V2). Consequently we have, by Lemma 3.1, g(#, @) 
=0 for every x1, x2, Vi, and V2; this proves the rest of the lemma, 
since U, Vi) — M3 (xe, V2) = M3(x1—xe, Vi- V2). 

To complete the proof of the multiplier rule we shall now prove 
that there is a linear continuous operation u=P(u) on MM to U such 
that ge(u)=J|[P(u)] for every uw. To this end define P(u) to be 
,[0., u, N(0., u)] and Q(u) to be #,[0,, u, N(0z, u)]. Therefore, by 
Corollary 2.3, 1=M;[0., P(u), Q(u)]; as a consequence it is evident 
that go(u) = ge(Ms[0., P(u), Q(u) which, by Lemma 3.2, is J[P(u)]_ 


THEOREM 2. There exist constants Io, C1, - - - , Cq, not all zero, and 
a linear, continuous, and real-valued functional J on UU such that 


(1) Lod, (Zo, Ho; + cud, Ho; — Nz, = 0 
for every p in M, and 
(2) lod sf (Zo, x) + Hu(Zo, Ho; x) — J(1[x, 0,, N(x, 0,)]) = 0 
for every xin 

From the previous discussion it is evident that 


ge(u) = lod f(Zo, Ho; + cxde(Zo, Ho; = J[P(u)] 


and that P(u) is ,[0., u, N(0., u) |. Equation (2), the “transversality 
condition,” follows from equation (3.1). For, 


go|M3(x, 0., 0.)] = J(P[M3(x, 0., 0.)]). 


Then making use, on the one hand, of the definitions of P and 9, and 
on the other hand of Corollary 2.3, we obtain the desired conclusion. 


4. Instances of the theory. We shall illustrate the generality of 
our theory by working out two special cases, the minima of functions 
of several variables and the problem of Bolza in the calculus of varia- 
tions. However the theory also includes the situations treated by 
Hahn, Graves, and the author, as well as many others. 

Consider first the problem of minimizing a real-valued function 
g(x1, - - -,Xm) of mreal variables in the class of points x9 = (x1, - - -,%m) 
satisfying equations h,(xo) =0, (k=1, ---, q). It is assumed that g 
has a minimum at a point Zo in this class and that g and h; have 
first differentials at this point. Then to indicate the applicability 
of our theory to this case we choose the sets Mt, MN, Meo, Mo, U, 
and % to be equal to the set of real numbers, ¥ to be the m-di- 
mensional euclidean space, and Xp to be the region of definition of 
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the functions g and fy, - - - , 4g. We then make the following defini- 
tions: Mo) =Ax(xo), (R=1, - - - , 9), f(X0, Mo) (x0), (Xo, wo) =0, 
¢1(X0, Mo, Yo) =o, and fo=0. Then it is clear that f(x9, uo) is a mini- 
mum at (Zo, Zo) in the class of all points (xo, wo) in the region which 
satisfy 


wo) = O, k=1,---,q, 
Ho = 0, 


and all previous assumptions are satisfied. 
Our multiplier rule then implies the existence of constants 
Io, C1, * - * » €g, not all zero, and a constant e such that 


log + = 0, edu = 0 
for every dx, ---, dxm and every dy; therefore we have e=0 and 
log Zo) + crhz(Zo) 0, i= 


which is the familiar form of the multiplier rule for this problem. 
The problem of minimizing a function f on a region Xo of a normed 
linear space in the class of points satisfying functional equations 
¥i(xo) =0, (R=1,---, q), can be treated in precisely the manner 
indicated above. 

Let us turn our attention to the problem of Bolza in the calculus 
of variations, where, for simplicity, we shall suppose that the admissi- 
ble arcs are of class C’. Here we are given an arc 9:(x), --- , 9,(x), 
<x <%.), which minimizes a functional 


I[x1, x,y, = yan), y(x2)] + f h(x, y(x), 


in the class of admissible arcs which satisfy certain side conditions 
galx, (x), a=1,---,r<s, 

and end conditions 
y(x1), x2, y(x2)] = 0, k= i,---,¢. 


For our purposes it is sufficient to restrict attention to the 
class Jo of all admissible sets of functions (yi(x),---, y,(x)), 
(%:—e <x <%,+¢e), where e is some fixed positive constant. The class 
which is composed of all sets (yi (x), - - - , y/ (x)) of derivatives of an 
element vo=(yi(x), - - , ¥s(x)) is Then if the transformation 


(4.1) yx) = + f (Oat 
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is performed upon the functionals J and ¥;, they become functionals 
(xo, wo) and Wi(xo, Of xo= (x1, Xe, y(%1)), Ho=(y’). Thus Xo is the 
set of all (x1, x2, yo, - Yos), Where | <e, | #—x2| <e, and 
(yor, - - * » ¥os) is the value of an admissible set at x = %. 

To complete our transformation of the problem of Bolza we define 
¢1(Xo, fo, Yo) to be the set of m functions 


¢1[x, vo(x), po(x) J, vo(x), uo(x) | 


and n(xo, Mo) to be the set of s functions of (xo, uo) which appear in 
the right members of equations (4.1). It is then a simple matter 
to verify assumptions A through D of §1, if it is assumed that 
the matrix ¢.,, of partial derivatives has rank m along the minimiz- 
ing arc. Assumption E can be verified without difficulty by making 
use of the theory of Volterra integral equations. 

The functional J appearing in our multiplier rule can be shown* to 
be representable as 


Now, if we define F to be Joh+)ada, identity (1) of Theorem 2 
yields the following equations: 


F,(x)dx + + + Wiis) = 0, 
(— 1)*F(%s) + + + | = 0, 
b = 1, 2; 


moreover identity (2) becomes 


By very obvious manipulations these equations can be put into the 
more familiar form 


— + + lodg + = 0, 
Fy -f F,dx + d;, 
which is the multiplier rule for the problem of Bolza. 


THE UNIVERSITY OF CHICAGO 


* See Hahn, loc. cit., pp. 544-546. 
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ON THE DERIVATIVE OF f(x)* 
H. S. WALL 


Let 41, ye, ¥3, - be defined recursively as follows: is the loga- 
rithmic derivative of a function y=f(x), and y,=D.y,.1, (v=2, 3, 
4,---). Then the successive derivatives y’, y’’, y’’’, - - - of y with re- 
spect to x are polynomials in y and the y,. In fact, y’=yy1, y’’ 
9” and 


(n) (n) vn 
where A®)...,. is a positive integer and the summation is taken for 
all non-negative integral solutions 1, v2, v3, -- - , ¥n of the equation 
(2) + + 3x3 +m, = 1. 


This statement may readily be proved by mathematical induction. 
The principal object of the present note is to prove the following 
theorem: 


THEOREM. The integer A =i om tn (1) is equal to the number of ways 
that n different objects can be placed in compartments, one in each of 
vy, compartments, two in each of ve compartments, three in each of vs 
compartments, --- , without regard to the order of arrangement of the 
compartments. 


1. Generalized binomial coefficients. Let k, m, n, (kn <m), be posi- 
tive integers, and denote by C,, the number of ways that kn objects 
can be selected from m objects and placed in compartments, & in 
each compartment, where no account is taken of the order of ar- 
rangement of the compartments. Thus C,”) is the binomial coefficient 
Cn,n=m!/[n\(m — n)!]. We have 


(k) 


or 


(k) 


(3) Cx, = m!/[n\(m — kn)\(k!)*]. 


This has meaning if m2kn. For special 0 values of the indices we 
shall consider C,®, to be defined by (3) by taking 0!=1. Thus if 
k=0, m=0, we have C{}=1. 


* Presented to the Society, September 5, 1936. 


= 


396 H. S. WALL [June 


If (2) holds, it will be seen that 


(2) (n) n! 


where ki} =n, ke ks We are to prove that 
this is the value of a=: -», in (1). For the proof we need the following 
identities which will be seen to hold for all values of m, n, k for which 
the symbols involved have been defined: 


(k+1) (k+1) (k+1) 
(k) 
(6) (n + 1)- atl = n 
Let it be remarked in passing that if Pr..«.=C%+C®%x+C®x? 
+ ---, then from (5) it follows that P,,4=Pr—1,e.+Cr1,rixPn—k,k. 


Also Pi sP 1: 


2. Derivation of the formula for Any: ..»,. Denote the sum in (1) 
by S,, and write S, as a polynomial in y,, S,= SO y,", where 
S © is independent of y;. We begin by showing that 
(7) Siz = 
We use induction on the subscript difference 1—j =k. From the rela- 
tion yS,=D, [ySn-1] it follows that 

. (1) (1) (1) 

(8) = y—1 + (1 + v)) y+1 + D, 


with the agreement that S;,;=0 if 7<0 or j7>7. Assuming that (7) 
holds for k<q we shall prove that it holds for k=g. Accordingly, 
we choose n, v, (0Sv<n), in (8) so that n»—v=g. Then, by our as- 
sumption, (8) may be written in the form 


(1) (1) (1) 


( S, (1 + v)¥2C -2,0 + 1,05 


(1) 
itg > 1s 
(9) Sn» } ay 


Replace 1 by n—», and v by 0 in (8), and eliminate D,S@,_,9 in 
(9). The result, by (6), is 


( (1) 


1 
(10) Da, = Soy. v—1 + ly 
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Hence 


i=0 i=0 


or, by (5) with k=0, S®?=C™S,_,.0, as was to be proved. 
We next put 


[m/p] 


= p = 2,3,4,--- 
Then the formulas 


hold for p=1. Assuming that they hold for p<k, k>1, we may then 
prove them for p=k. To do this, put y=0 and p=k—1 in (12), and 
equate coefficients of like powers of y;. The result is the equation (12) 
with p=k. Thus (12) is true when p=k, and in particular Ss 
= Hence by (5) we find that so that 
(11) holds for =k provided m—kv=0. The proof of (11) for p=k 
may now be carried out along the lines of the proof of (7). with in- 
duction, in this case, on the difference m— kp. 
After (11) has been proved it follows at once by (4) that 


(n) n! 
(13) ‘ 


3. Application. In conclusion I shall give examples to illustrate the 
application of the foregoing result. 


EXAMPLE 1. Let @,,, denote the number of ways that m different 
objects can be distributed among k compartments, where no account 
is taken of the order of arrangement of the compartments, and at 
least one object is placed in each compartment. Then elementary 
considerations will show that a,..=Rdn_1,.4+@n-1,.4-1. Put 


y= eter = L,(t)x’/v! 
v=0 
Then by (1) we find that L,(t) -- nl”. 
Put It follows that gi(y)=1/(1—y), 
gx(9) = (k=2, 3,4, - - - ), or ge(y) =1/[(1—y) (1 —2y) 
--+ (1—ky)]. Hence L,(1), the number of ways that » different ob- 
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jects can be distributed among z or fewer compartments, is the coeffi- 
cient of y” in the power series P(y) for the function 
y y* 
where m=n. The number L,(1) is also the coefficient of x"/n! in the 


power series for the function e“*-. 
The a,,, are given explicitly by the formula 
(- 1 k+1 
EXAMPLE 2. Put y=(1+<)-‘ in (1) and then set x=0. There re- 
sults this identity: 
+ 1)¢ + 2)---¢+n- 1) 


n! 


(14) 


priteet 


(1-2- -v)(2-4- - 2ve) (n-2n- - Nn) 


where the summation is taken as in (1). On putting ¢=1 in (14) we 
obtain the following theorem :* 


THEOREM. Form a partition 2 n by gop at most one a from 
each of the progressions 1, 2, 3, 2, Oy 
Multiply together the terms of to including 
the integer chosen. Let the products so formed be a, b, c,---. Then 
[1/(a-b-c- --- )]=1, where the sum is taken for all such 
of n. 


EXAMPLE 3. If we differentiate the members of (14) with respect to 
t and then set ¢=1, we get the formula 
1 1 1 
2 3 n 


> (vy + ve +) 


This may likewise be interpreted as a theorem on partitions of n. 


NORTHWESTERN UNIVERSITY 


* facobi, Zur combinatorischen Analysis, Crelle’s Journal, vol. 22 (1841), pp. 372- 
374. 
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THE GEOMETRY OF THE WHIRL-MOTION GROUP G;: 
ELEMENTARY INVARIANTS* 


EDWARD KASNER AND JOHN DE CICCO 


In this paper we shall study the elementary geometry of the ori- 
ented lineal elements of the plane with respect to the whirl-motion 
group Gs. We give results in addition to those found in a paper by 
Kasner, The group of turns and slides and the geometry of turbines, 
published in 1911 in the American Journal of Mathematics (vol. 33, 
pp. 193-202) and the paper by the authors, The geometry of turbines, 
flat fields, and differential equations, published in 1937 in the Ameri- 
can Journal of Mathematics (vol. 59, pp. 545-563). The present paper 
can be read independently of the two earlier papers and contains the 
foundation of a new geometry of differential elements. 

We begin by considering certain simple operations or transforma- 
tions on the oriented lineal elements of the plane. A turn T, converts 
each element into one having the same point and a direction making 
a fixed angle a with the original direction. By a slide S; the line of 
the element remains the same and the point moves along the line a 
fixed distance k. These transformations together generate a continu- 
ous group of three parameters which we call the group of whirl trans- 
formations. It is noted that the group of whirls is isomorphic to the 
group of motions. These two three-parameter groups are commuta- 
tive and together generate a new group of six parameters which we 
term the whirl-motion group Gs. It is our purpose to study the simple 
geometry of this group Gs. We find the fundamental invariants. 

We define ~'! elements to be a series of elements; this includes a 
union or curve as a special case. A set of ©? elements is said to forma 
field of elements, which corresponds to a differential equation of the 
first order F(x, y, y’) =0. A turbine is the series which is obtained by 
applying a turn 7, to the elements of an oriented circle (which may 
be a null circle). A flat field is the field that is obtained from the total- 
ity of all elements which are determined by the set of all oriented 
circles containing a given element. We obtain the elementary rela- 
tionships between elements, turbines, and flat fields. 

For the analytic representation, it will be convenient to define an 
element by the coordinates (u, v, w) where v is the perpendicular from 
the origin, u is the angle between the perpendicular and the initial 


* Presented to the Society, December 30, 1937. 
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line, and w is the distance between the foot of the perpendicular and 
the point of the element. 


1. The whirl-motion group G, of transformations.* The equations 
of the turn T, are 


(1) #=u+a,i =vceosa+wsina, = —vsina+weosa. 
The equations of the slide S; are 
(2) i=, 2, w=wtek. 
Since any whirl may be expressed in the form 7;S;T., the equations 
of any whirl are 
+ ksinB, 
@ = — vsin (a+ 8B) + woos (a+ + k cos 8B. 


The only contact transformations of the group of whirl transforma- 
tions are 


(3) 


T_+/2S aT T aT 


The first represents a dilatation Dg and the second, which may be 
written 7,Da, represents a dilatation accompanied by a reversal of 
orientation. 

It is then seen that any whirl transformation may be given by 


SiDaT a; 


so that the equations of any whirl transformation may be put in the 
form 


=ut+a, 


(4) 


 k. 


It is observed that the group of whirl transformations and the 
group of motions are isomorphic. These two groups are commutative 
and together generate a continuous group of six parameters which 
we call the whirl-motion group Gg. 

The equations of any motion are 


a=utn, v=v+ypcosu+vsinu, 
w= w—psinu+rcosu. 


From (4) and (5), we find that the equations of any whirl-motion 
transformation are 


* See Kasner, The group of turns and slides and the geometry of turbines, loc. cit. 
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(6) 
= —vsina+ weosa 


This is our fundamental six-parameter group Ge, containing the whirl 
group W3;, and the motion group M; as three-parameter subgroups. 


2. The invariants for two elements. The angle between two ele- 
ments is the angle between their lines. The distance between two 
parallel elements is the distance between their points. 


THEOREM 1. Under Ge we find that (1) two non-parallel elements 
possess the angle between them as the unique invariant, and that (2) 
parallel elements are transformed into parallel elements in such a way 
that the distance between them 1s preserved. 


3. The turbine. A turbine is the series which is obtained by apply- 
ing a turn 7, to the elements of an oriented circle. 
The equations of a turbine are 


v=acosu+bdsinu+r, w= —asinu+bcosu+s. 


We call (a, b, r, s) a set of turbine coordinates. 

Under Ge, any turbine is converted into a turbine. Thus Gg in- 
duces a correspondence between the turbines of the plane. By (6) we 
find that this correspondence is given by the equations 


@=acosdrX — bsink +ywcosdA — vsind, 
+ pcos kX, 
7) 7 =rcosa+ssina+d, 

$=-rsina+scosat+k. 


From the equations of a turbine, we obtain the following theorem: 
THEOREM 2. Two non-parallel elements determine a unique turbine. 


THEOREM 3. The necessary and sufficient condition that two distinct 
turbines possess a common element (in which case they have only one 
common element) is 


(a2 — a)? + (b2 — = (re — 13)? + (Se — 51)?. 
THEOREM 4. Two turbines have two and only two invariants, namely, 
(a2 = a)? + (b2 — b;)? and (re — 17)? + (se — $;)?. 


From this and Theorem 3, it follows that two intersecting turbines 
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possess the squared distance between their centers as the unique in- 
variant. 


THEOREM 5. A given element and a given turbine possess the unique 
invariant 


ja—(v—r) cos u+(w—s) sin u]?+[b—(v—7) sin u—(w—s) cos u]?. 


From this it is seen that a given element E and a given turbine T 
possess as the unique invariant the squared distance between the ele- 
ment E and that element on T which is parallel to E. 


4. Conjugate turbines. Two turbines T and T are said to be con- 
jugate if they have the same circle as point locus and the elements of 
the two turbines are symmetrically related to the elements of the 
circle. 

The two turbines T(d, 5, #, 5) and T(a, b, r, s) are conjugate if and 
only if = 

ag=a, 

THEOREM 6. The conjugate turbines of two given turbines do or do not 
possess a common element according as the two given turbines do or do 
not possess a common element. 


5. The flat field. The totality of elements determined by the set 
of all circles which contain a given element is called a flat field. The 
given element is called the center (or central element) of the flat field. 

The equation of a flat field is 


w= — (v — d) cot 3a — u) — ®, 


where (i, 3, @) are the hessian coordinates of the central element. 

Under the whirl-motion group Gg, any flat field is converted into a 
flat field. Thus G, induces a transformation between the flat fields of 
the plane. By equation (6), this induced correspondence is given by 
the equations 


ai=ti—a+n, 
(8) idcosa — sin a+ucos(é — a) sin — a) +d, 
= dsina+ cos a — psin (4 — a) + cos — a) — k. 


It is seen that (8) is a whirl-motion transformation. We call it the 
conjugate whirl-motion transformation WM of the whirl-motion trans- 
formation WM given by equation (6). 


THEOREM 7. Three elements which are not all on the same turbine and 
which do not all possess the same direction determine a unique flat field. 
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This follows from the equation of a flat field given above. 


THEOREM 8. The turbines which are contained in a flat field are those 
whose conjugate turbines contain the central element of the flat field. 


THEOREM 9. Two flat fields whose central elements are not parallel 
always have a common turbine. Three flat fields, no two of whose central 
elements are parallel, have a common element or else a common turbine. 


DEFINITIONS. The angle between two given flat fields is the angle 
between their central elements. Two flat fields are said to be parallel 
or supplementary (anti-parallel) according as the angle between them 
is 0 or x. The distance between two parallel flat fields is defined to be 
the distance between their central elements. 


THEOREM 10. Under Ge, we find that (1) two non-parallel flat fields 
possess the angle between them as the unique invariant, and that (2) two 
parallel flat fields are transformed into parallel flat fields in such a way 
that the distance between them is preserved. 


THEOREM 11. A given element and a given flat field possess the unique 
invariant 


| @ - V) cos ere (w + W) sin 5 i. 


From this, we find the following geometric interpretation of the in- 
variant: Let E be the given element and F the given flat field. Let the 
central element of F be denoted by G. On the oriented line of E, 
construct the element G which is in the flat field F. Let / be the line 
connecting the points of G and G. The square of the perpendicular 
distance from the point of E to the line / is our invariant. 


THEOREM 12. A given turbine and a given flat field possess the unique 
invariant 


[a—(d—r) cos i+(@+s) sin [b—(5—1) sin —(H+s) cos 


From this we find the following geometric interpretation of our in- 
variant. Let T be the given turbine and F the given flat field, let T 
be the conjugate turbine of T, and let G be the central element of F. 
Then our invariant is the square of the distance from the element G 
to the turbine T. (See construction in Theorem 5.) 
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A NOTE ON NUMBERS OF THE FORM 


ROBERT M. THRALL 


It has long been known that the only sets of values of a, 8 which 
give rise to universal forms a?+ab?+c?+afd? are a=1, B=1, 2, 3; 
a=2,8=2, 3,4, 5. We give here a theorem from which the universal 
character of all these multiplicative universal forms can be readily es- 
tablished. The methods and arguments used are completely algebraic 
in character. This note is closely allied to R. D. Carmichael’s paper, 
Proof that every positive integer is a sum of four integral squares.* We 
use the formulas on the first page of his proof without recording them 
here. 

We remark that every prime number ? is a divisor of the form 
a?+ab?+8c?+afd? in which a, b, c, d are relatively prime. For if p 
is prime to a, a theorem from the theory of quadratic residues states 
that p isadivisor of a?+ab?+ 8. If p isnot primeto a, takea =c=d=0, 
b=1. In both cases the bases are seen to be relatively prime. 

We now require that a, B be positive integers with B2a>0. The 
results obtainable also hold for a=0, which case has been treated by 
R. D. Carmichael. 


THEOREM 1. For a<3 and every prime number p>B there exists a 
positive integer q < [48/(3—a) such that p 


For p<5 we verify the truth of the theorem directly. Henceforth 
consider p=5. Since p>S 2a, from the above remark we have a q’ 
such that pg’ =a?+ab?+ 8, where a and b may evidently be taken 
less than p/2. Hence, pg’ S(1+a)p?/4+6 or g’S(1+a)p/4+8/p 
<3p/4+1; that is, g’< pif p25, which is the case now under con- 
sideration. 

Let g be the smallest positive integer such that pg=a?+ab?+ 8c? 
+aBd?. From the above we have g<p. Also a, 5, c, d are relatively 
prime, for otherwise the square of their greatest common divisor 
would divide g, leaving pq: in the same form as pg with qi <q contrary 
to the hypothesis that g has the smallest value possible. 

We now consider the forms in the parentheses of the formula that 
would correspond to Carmichael’s (2): 


* Duke Mathematical Journal, vol. 2 (1936), pp. 243-245. 
¢ American Mathematical Monthly, vol. 44 (1937), pp. 81-86. 
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A; = ax — aby — Bcz — afdt, 
Ag = bx + ay + Bdz — Bet," 
A; = cx — ady + az + abt, 
A,g=dx+cy — bz + at. 


The expression A, is invariant under any change of x, y, 2, ¢ of the 


form: 
x + asc + + a30, 


x 
y’ 7 = asad, 
t 


2+ acd + asc + aca, 
t — azsd — asc + acd, 


where the a’s are arbitrary. Under this transformation A; becomes 
A3’ = A3 + + ad?) + (a2 — aas)(cb + ad) 
+ (a3 + as)(ac — abd) + ag(a? + abd?). 


We may choose x =¢=0, y=b, z=c so that A,=0. If a prime p di- 
vides a number of the form a?+ab?, (a, b relatively prime, a=1 or 2), 
then ? is itself of that form, hence g=1 < [48/(3—a) ]*/?. For other 
primes the greatest common divisor 6 of a?+ab? and c*+ad? must be 
1. For if has a prime power factor then (a?+ab*)/p," =a? +ab? ; 
(c?+ad?)/p," =c?+ad? (this follows from the theory of divisors of 
the forms a?+8? and a?+25?); and (g/f,") has the form considered. 
We can therefore choose a, a, (taking a2=a;=a,=a;=0) so that 
A;'=1. Then in (2)* take ¢=p=0 and X, yw such that the bases in the 
first two braces are numerically less than or equal to g/2. Now sup- 
pose that g > [48/(3—a) ]"”?; that is, that B<q?(3—a)/4. The form A 
is replaced by A’ with the new choice of x’, y’, 2’, t’, and 


<A’ 2 = 


Hence 0 <A’ <q; but in view of Carmichael’s formula (6) this contra- 
dicts the minimal nature of g. Therefore, gS [48/(3—a) ]"?, and the 
theorem is proved. 

We now give a corollary which is useful in some applications of the 
theorem. 


Coro.iary 1. If and a+B+1<%9, then gS2. 


For, referring to (2) with g =3, we see that by choosing x’, y’, 2’, z’, 


* References are to equations in Carmichael’s proof. 
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ut, A, p, o as in the proof of Theorem 1 we have 3A’<1+a+8<9. 
Therefore A’ <3. But A’ is an integer, and if it is less than three this 
contradicts the minimal character of gq. 

A form ax?+by?+cz?+d?? is called universal when it represents all 
positive integers. From the multiplicative property of the form 
a*?+ab?+8c?+ afd? we see that if it represents every prime number it 
is universal. If the number gq, defined in the preceding theorem, is 
one, then the form does represent every prime and therefore every 
positive integer. We turn to the consideration of particular values for 
a, 8. In the following paragraphs the theory of the forms a?+ 8c? 
will be assumed. 

The case a=6 =1 gives g=2'/?. Hence g=1, and we have the fol- 
lowing theorem: 


THEOREM 2. Every positive integer is a sum of four integral squares. 


That the forms a?+ab?+6c?+afd?, with (a, B)=(1, 2), (1, 3), 
(2, 2), (2, 3), (2, 4) are universal follows readily by application of 
Theorem 1 and its corollary. We shall treat just one typical case, 
a=2, 8=3, among these. 

Ifa=2,8 =3, we have S (12)! org $3. But, sincea+B+1=6<9, 
the corollary applies and g <2. If 2b=a?+26?+3c?+ 6d’, then a and c 
are of like parity. If they are both even, p=b?+2(a/2)?+3d?+6(c/2)?. 
If they are both odd, by proper choice of sign for a and c we have 
where A+3u =a, and p=b?+2\?+ 3d? 
+6y?. We have now proved the following theorem: 


THEOREM 3. Every positive integer is representable in the form 
a?+2b?+ 3c?+ 6d?. 
The case a=2, B=5 gives gS(20)'” or gS4. But if 4p=a?+20? 
+5c?+10d?, (p an odd prime), then a and c must both be even, and 


2p=b?+2(a/2)?+5d?+10(c/2)?. Hence But, since a+6+1 
=8 <9, the corollary gives g <2, and 


(1) 2p = a? + 26? + 5c? + 10d?, 


where a and c must be of like parity. If they are even, p=b?+2(a/2)? 
+5d?+10(c/2)?. Hence suppose 


(2) a,c odd. 
Then by taking residues (mod 4) in (1), 
(3) b = d (mod 2). 


Now by (1) 
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a+ 2b + 5c + 10d\? a—b+ 5c — 5d\? 
r=( 6 +2/ 6 ) 


2b — c — 2d\? —~b—c+d\? 


(4)* 


In view of (2) and (3), the numerators in (4) are all even. Then, 
unless exactly three of a, b, c, d are divisible by 3, we can choose 
signs for a, b, c, d so that 


(5) a—b—c+d=0 (mod 3). 


Then all the other numerators in (4) are divisible by 3. 
In the exceptional case either a and b or c and d are divisible by 3. 
But the identity 


(6) 9(A? + 2B?) = (A + 4B)? + 2(2A F B)? 


(repeated if necessary) shows that any multiple of 3 of the form 
x?-+2y? can be expressed in that form with x, y prime to 3. Then (5) 
can be verified as above, and g = 1. We have now proved the following 
theorem: 


THEOREM 4. Every positive integer is representable in the form 
a?+ 2b? + 5c?+ 10d?. 


UNIVERSITY OF ILLINOIS 


A MOMENT-GENERATING FUNCTION WHICH IS USEFUL 
IN SOLVING CERTAIN MATCHING PROBLEMS 


EDWIN G. OLDS 


1. Introduction. In a book published several years ago, Fryf de- 
voted considerable attention to various aspects of a problem which 
he called, “the psychic research probiem.” His introductory problem 
is the following: “A spiritualistic medium claims to be able to tell the 


* Formula (4) and the rest of the proof of this theorem were suggested by 
Gordon Pall. 

Tt Presented to the Society and The Institute of Mathematical Statistics, Decem- 
ber 30, 1937. 

t T. C. Fry, Probability and Its Engineering Uses, Van Nostrand, New York, 1928, 
pp. 41-77. 
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color of a playing card without seeing it. In order to test her claims 
an experiment is conducted with four red and four black cards. These 
cards are thoroughly shuffled and placed face down on the table. The 
medium is told that there are four red and four black cards, but pre- 
sumably knows nothing as to their arrangement. The experimenter 
picks up a card and, without either looking at it himself or showing 
it to the medium, asks its color. If she answers ‘red,’ he places it at 
one side of the table. If she answers ‘black,’ he places it on the other 
side of the table. This process is repeated until all cards are exhausted. 
If the medium does not have the ability which she claims to possess, 
what is the chance that there will be just one black card in the pile 
that should be red?” 

Fry solves this problem by considering the number of chance orders 
of red and black which would match the order in which the medium 
calls the cards to the extent of producing exactly three red pairs. He 
also proposes and solves various other problems of this general type. 
A problem like Fry’s, but concerned with tea-tasting, is thoroughly 
discussed by Fisher.* 

Chapmanf in 1931, contributed the solutions of several problems 
related to the expected number of correct matchings if a deck of ¢ 
different cards is placed in correspondence with a given target deck 
of the same kind. One of his results is the function 


which gives the probability of s correct matchings. 

A more general problem has been considered by Greenwood.{ He 
considers matching a target deck containing ¢ sets of identical cards, 
with s; cards in the ith set, ({=1, 2, - - - , t). For the case where the 
number of cards is the same for each set, he obtains the mean and 
variance of the distribution of the numbers of correct matchings as s 
and s*(¢—1)/(st—1), respectively. The purpose of the present paper 
is to develop a moment-generating function and derive the first four 
moments for the theoretical distribution of numbers of correct match- 


*R. A. Fisher, The Design of Experiments, Oliver and Boyd, Edinburgh, 1935, 
chap. 2. 

¢ D. W. Chapman, The statistics of the method of correct matchings, American 
Journal of Psychology, vol. 46 (1934), pp. 287-298. 

t J. A. Greenwood, Variance of a general matching problem, Annals of Mathe- 
matical Statistics, vol. 9 (1938), pp. 56-59. See also M. S. Bartlett, Properties of 
sufficiency and statistical tests, Proceedings of the Royal Society, (A), vol. 160, pp. 
268-282 (note especially the discussion of the / Xm contingency table, pp. 278-280). 
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ings for this case where the number of cards is the same for each set.* 


2. Discussion of Fry’s problem. Before taking up the general case 
of ¢ sets of s cards each, let us consider a solution of Fry’s problem, 
which we quoted above. 

Suppose the cards are tagged Rj, Re, Rs, Rs, Bi, Be, Bs, By. Let us 
consider the order in which the medium calls the cards as constituting 
the order in which a target deck is laid out. When the deck used for 
matching is laid out in correspondence to the target deck, any of 
the cards of the matching deck may fall below any one of the cards of 
the target deck. This is illustrated in the array below, where the col- 
umns correspond to the cards of the target deck and the entries in 


R,| Re| Rs| Bs | Be| Bs| Be R,| Re| Rs| Rs| Bi | Bs | Bs 
Re| Re| Re| Re| Re| Re| Re 
Rs| Rs| Rs| Rs| Rs| Rs| Rs| Rs 
Ra| Ra| Ral Re 1 
B,| B,| By| By| By| By 
B,| Bz| Bz| Bz| By| Be 
B;| Bs| Bs| Bs| Bs| Bs| Bs| Ba| 
By| Bs| Bs| Bs| Bs| Bs| Bs| Be 
Array I Array II 


the columns correspond to the cards of the matching deck. Now, 
whenever R; appears in a column labeled Rj, Re, Rs, or Ry, this repre- 
sents a matched pair; and likewise other reds in the red columns and 
the blacks in the black columns represent matched pairs. In every 
such case let us replace the symbol by x, and in all other cases by 1. 
In this way we form the second array. 

Now, if we make all possible choices of one element from each 
column of the first array, with no two from any row, we have all of 
the 8! orders in which the cards of the matching deck may appear. 
Let us suppose that one of these is R,R2B,B;R;R,B2B,. This corre- 
sponds to x, x, 1, 1, 1, 1, x, x in the second array, or, if we drop the 


* T. E. Sterne, The solution of a problem in probability, Science, vol. 86, pp. 500-— 
501, has calculated the first four moments for the case s =5, t=5; also note numerical 
results given by E. V. Huntington, Exact probabilities in certain card-matching prob- 
lems, Science, vol. 86, pp. 499-500; and results given by C. E. Kellogg, New evidence(?) 
for extra-sensory perception, Scientific Monthly, vol. 45 (1937), pp. 331-341. 
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commas and multiply, it corresponds to x‘, where the exponent indi- 
cates that we have exactly four matched pairs. Forming every possi- 
ble product from the second array, using only one element from each 
column and no two from any row, we have 8! terms of various de- 
grees. Collecting like powers we have 


= (41228 + 16(4!)2x* + 36(4!)2x4 + 16(4!)22? + (412, 


which means that, of all the possible arrangements, (4!)? produce 
eight matched pairs, 16(4!)? produce six matched pairs, for example. 
If we divide the coefficients by 8! we have the probabilities. Thus 
Fry’s problem, which requires the probability for six matched pairs, 
has the answer 16(4!)?/8! or 8/35. 

If we denote by f(r) the frequency of r matched pairs, we may write 
the above polynomial in the form 


¢(x) = f(r)x". 


Then, obviously, ¢(x)/8! is a generating function for factorial mo- 
ments. Denoting the factorial moments by K, where 


K;=— if =r(r—1)---(r —i+1), 
8! 
we have 
Ko = ¢(1)/8!, Ki = ¢’(1)/8!, Ki; = ¢*(1)/8!. 


The form of the function ¢(x) needs detailed consideration. This 
polynomial was obtained from the array of x’s and 1’s by adding all 
the terms formed by choosing one element from each column but no 
two from any row. Such an array has been called a permanent,* and 
the symbol +||+ has been used for it. The properties of the perma- 
nent which will be useful in this paper are given in the next section. 


3. The permanent. It is obvious that the permanent has properties 
analogous to those properties of the determinant which do not depend 
on the fact that the determinant vanishes when two rows (or two 
columns) are identical. For example, the permanent can be expanded 
in terms of minor permanents. Of particular importance is the fact 
that the derivative of a permanent in one variable is equal to the sum 


* T. Muir, revised by W. H. Metzler, A Treatise on the Theory of Determinants, 
Muir, Albany, 1930, p. 19. 
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of the permanents formed by replacing, one column at a time, the 
elements of each column by their derivatives. 


4. Application of the permanent. As noted above, we may identify 
(x) as a permanent, namely 


behit filet 


Obviously, ¢(1) =8!. The first term of $’(x) is 


ot 


and it is easy to see that each of the other seven terms is equivalent 
to this one. Therefore 


= 4-8-}1 11 and ¢(1) = 4-8-7! =4-8!. 


In the permanent above, one of the blocks of x’s has one less col- 
umn and one less row than before. It is convenient, then, to denote 
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this permanent by [1]. Using the same system, we will denote by 
[21] the permanent, where one block of x’s is reduced by two rows 
and two columns and the other by one row and one column. 

It can be readily verified that ¢’’(x) =4-8(3-3[2]+4-4[11]) and 
=4-8(3?+4?*) -6!=(4-5?-8!)/7. The succeeding derivatives may 
be written in similar fashion. 

Using these results, we may compute the factorial moments and, 
from them, the moments about the mean. 


5. The solution of the general problem. For the general problem 
we have ¢ groups (instead of 2) of s like cards each (instead of 4). 
The general method of obtaining the moment-generating function, 
and from it the moments, is the same as was used for the special 
case, except that the work becomes somewhat more intricate and 
arduous. In the general case the function ¢(x) is a permanent of n 
rows and z columns, (st=m), and has ¢ blocks of s x’s each. Then 

¢(x) = [0], 
¢'(x) = 
(x) = (s)%[2] + s# [17], 
= (s))%4[3] + [21] + [13], 
= (s)2¢[4] 4+ [31] + 3(s(2))4¢(2) [22] 
+ 6(s?)s2)%4@ [212] + [14], 


and, in general, 


where the notation is explained as follows: 
1. ~;7* indicates that p; appears 7; times in the partition of m. 
2. The summation is to include all possible partitions of m. 
3. For each term, 
m! 
C.= 


The work of Dwyer,* in a recent paper, proved to be of consider- 
able assistance in obtaining a compact expression for the general de- 
rivative, as given above. It should be noted also that the numerical 
coefficients of the partitions which occur in our derivatives are given 


* P.S. Dwyer, Moments of any rational integral isobaric sample moment function, 
Annals of Mathematical Statistics, vol. 8, pp. 21-65. 
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in his Table II.* For example, for the fourth derivative we might have 
used his “weight =4” and read the proper numerical coefficients from 
the first column, namely, 1, 4, 3, 6, 1. 

To evaluate these derivatives of ¢(x) at 1 =1 it is sufficient to no- 
tice that, if the sum of the bracketed numbers is j, the permanent, 
represented by the bracket, has the value (m— j)! for x=1. For ex- 
ample, each of the bracketed quantities in ¢‘¥(x) becomes (n—4)! 
when x is set equal to 1. 

The rest of the work necessary to calculate the moments about the 
mean is quite direct, although tedious. We obtain the factorial mo- 
ments by dividing ¢‘(x) by m!; then, by means of well known rela- 
tions, we obtain the moments about the mean. Representing the 
moments about the mean by y; and replacing t by ”/s, we obtain the 
mean and higher moments as follows: 


s(n — s) 
s(n — s)(n — 2s) 
s(n — s) 
= {(n — 2s)(m — 3s)(3s + 1) 


+ (s — 1)(12ms — nm — 18s? — 6s)}. 


Comparisons with the articles cited above indicate complete agree- 
ment between the results given for various special cases and the val- 
ues obtained by the use of the formulas derived in this paper. 

In conclusion, it should be mentioned that this general method can 
be applied to the more comprehensive case where the numbers in the ¢ 
sets of the matching deck are not necessarily the same but may differ 
from set to set. Curiously, the mean is still s, but the variance is de- 
creased. 


CARNEGIE INSTITUTE OF TECHNOLOGY 


* Loc. cit., pp. 30-32. 


A THEOREM ON INTERIOR TRANSFORMATIONS* 
G. T. WHYBURN 


In an earlier paperf it has been shown that if 7(A) =B is a light 
interior transformation, that is, a continuous transformation map- 
ping open sets into open sets and mapping no continuum into a 
single point, and if A is compact, then for every simple arc pq in B 
and any poe7—'(p) there exists a simple arc pogo in A which maps 
topologically onto pg under T. 

In this paper the following extension to arbitrary dendrites (that 
is, locally connected continua containing no simple closed curves) in 
B will be made. 


THEOREM. Let T(A)=B be interior and light, where A is compact. 
For any dendrite D in B and any xoeT—'(D) there exists a dendrite E 
in A containing xo which maps topologically onto D under T. 


ProorF. Since the transformation T7—'(D) =D is interior,{ clearly 
there is no loss of generality in assuming that D=B. Now let H 
=> paqi be an arc-development of B; that is, for each 7>0, 
and (i=1, 2,---), and all points of H—H 
are end points of B, where we may suppose that po>=T7(x0). Now 
by the result quoted above, there exists an arc xoyo such that 
T(xovo) = Pogo and T is topological on xoyo. Likewise there exists an 
arc xyyi in A such that 7(x1y1)=figq, T is topological on x1, and 
furthermore so that y;=7~'(q:)-xoyo. Similarly there is an arc X2ye 
contained in A so that T(x2y2) = poge, T is topological on xeye, and ye 
=T-(q2)->-ixxy;. Continuing this process indefinitely we obtain 
K =>°xvy;, so that T(K) =H, and T is topological on K. Clearly K is 
a continuum and 7(K) =B. Hence our proof will be complete as soon 
as we show that T is 1 to 1 on K, or, what amounts to the same thing, 
that for each peB, T-'(p) -K reduces to a single point. 

Now there exists a monotone decreasing sequence of connected 
neighborhoods Vi, V2, V3, - - - of pin B with 6(V;)—0. Furthermore, 


* Presented to the Society, December 29, 1937, under the title Interior transforma- 
tions on certain curves. 

t See my paper in the Duke Mathematical Journal, vol. 3 (1937), p. 377, Theorem 
4.1. Compare with Stoilow, Annales Scientifiques de l’Ecole Normale Supérieure, 
vol. 63 (1928), pp. 347-382; and Montgomery, Transactions of this Society, vol. 42 
(1937), pp. 328-329. 

t See my paper, loc. cit., p. 370, Lemma 1.2. 


€ 
414 G. T. WHYBURN [June 


1938] INTERIOR TRANSFORMATIONS 415 


since H—H contains only end points of B, it follows that H- Vi, 
H- V2, - - - are connected sets which, of course, are open in H. Hence 
K-T“(Vi), K-T-(V2),- ++ is a monotone decreasing sequence of 
connected sets in K which are open in K. Then let L=lim K-T—(V;). 
Since V;>p, we have Lc T—(p). But, since L is necessarily con- 
nected and 7~'(p) is totally disconnected, L must reduce to a single 
point ge7—'(p). Hence -K =q, and our theorem is established. 

That a similar conclusion cannot be obtained, if we permit D to 
contain simple closed curves, is readily seen, since it fails to hold even 
in the simple transformation w =z? of the circle |z| =1 into the circle 
| w| =1. Still more striking, however, is the following: 


EXAMPLE. There exist two connected graphs A and B and an interior 
transformation T(A)=B, where A contains no subset homeomorphic 
with B. 


Let J =a,bificid2befece be a simple closed curve in a plane, where 
the points a;, b,, - - - are cyclically ordered on J as indicated. Let 
a;d:f2 and b,e,c, be disjoint arcs lying within J except for their end 
points, and let de, be an arc within J having only its end points in 
common with a:d,f2+,e,¢. Similarly, let fidea, and beeece be disjoint 
arcs lying, except for their end points, without J, and let d2¢. be an arc 
without J having just its end points in common with f,dea2+ bee2¢e. 
Finally let A be the graph thus constructed, that is, 


A =J + aydife + + dye: + fidede + + deez 
= + Difi + + + byes + + aids + dife + dies 
+ + + foce + Cod: + bree + + + + deez. 


Then A is a graph of eighteen edges as indicated in the latter sum. 
We now construct a graph B in 3-space as follows: Let @ be a 6- 
curve bac+bec+bfc in a plane 7. Let d be a point not in 7, aad let 
da, de, and df be arcs having just d in common and having just a, eé, 
and f, respectively, in common with @. 
If we now let 


B=0+da+de+df 
= ab + bf + fe + ca + be + ec + ad + df + de, 
then B is a graph of nine edges as here indicated. 
Now let 7(A)=B be the transformation which maps 4,1, debe 


topologically into ab, bifi, and bef. topologically into bf, and so on. In 
general T maps any edge x;y; of A topologically, preserving end 
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points, into the edge xy of B. Then T is an interior transformation. 
In fact T is a local homeomorphism which is 2 to 1. 

Since A is a planar graph, whereas B is non-planar (B is, in fact, one 
of the two well known Kuratowski primitive skew curves), clearly A 
contains no subset homeomorphic with B. Incidentally this example 
shows that planarity is not an interior property (that is, it is not in- 
variant under interior transformations). 


UNIVERSITY OF VIRGINIA 


ON THE TRANSFORMATION GROUP FOR DIABOLIC 
MAGIC SQUARES OF ORDER FOUR* 


BARKLEY ROSSER AND R. J. WALKER 


This paper concerns only magic squares of order four, and all state- 
ments of the paper are to be construed as applying only to magic 
squares of order four. 

One says that 


a b c d 

e f g h 
(1) 

1 j k l 

m n p 
is a diabolic (or pan-diagonal or Nasik) magic square if a, b,---,p 
are 1, 2, - - - , 16 in some order, and each row, column, and diagonal 


adds up to 34. This is to include broken diagonals such as 1, f, c, p, 
or c, h, i, n. A diabolic magic square clearly remains diabolic if sub- 
jected to the following transformations: 

A. Reflection about the a, f, k, p diagonal. 

B. Rotation through 90° counter-clockwise. 

C. Putting the first column last. 

D. Putting the first row last. 


For many purposes it is convenient to consider a diabolic magic 


* Presented to the Society, December 30, 1937. 
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square as mapped on a torus. That is, the square is bent and stretched 
until the top and bottom edges coincide and the right and left edges 
coincide. Then a diagonal like b, g, 1, m is continuous and not broken. 


THEOREM 1. The four elements of any square of order two of a diabolic 
magic square of order four mapped on a torus add up to 34. 


This includes squares like a, d, m, p or e, h, i, l. 

Proor. Add together the two main diagonals and the second and 
third rows and subtract the first and fourth columns. This gives 
f+g+j+k=34. By use of transformations C and D, this proof can be 
applied to any square of order two. 


THEOREM 2. If a diabolic magic square of order four be mapped on a 
torus, then the sum of any element and the element which is two distant 
from it along a diagonal is 17. 


Note that if one starts at any element and moves two elements 
along any diagonal, one always stops at the same element regardless 
of the direction of the diagonal. Two elements so situated will here be 
called antipodal. 


Proor. Add together the first and the third rows, then add the 
first and third columns. and the diagonals a, f, k, p and a, n, k, h; and 
subtract the diagonals e, b, 0,1; 7, f, c, p; m, j, g, d; and 1, n, c, h. This 
gives a+k =17. By use of transformations C and D, this proof can be 
applied to any pair of antipodal elements. 

By use of these two theorems, one can readily see that a diabolic 
magic square remains diabolic if subjected to the following transfor- 
mation: 


E. The transformation of (1) into 


a d h e 
b c g f 
n 0 k j 
| m p l 1 


THEOREM 3. All diabolic magic squares of order four can be derived 
from a single one by successive applications of A, B, C, D, and E. 
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ProoF. By use of C and D, 1 can be brought into the upper left- 
hand corner. Thus there is no loss of generality in taking a=1, k =16. 
Then c and d can at most be 14 and 15, so that c+d<29. Hence b2=4. 
Similarly c24, d24, e24, 124, m24. Since a+b+e+f=34, 
at+d+e+h=34, a+b+m+n=34, and a+d+m+p=34 by Theo- 
rem 1, one can similarly get f24, h2=4, n24, p24. Thus 2 and 3 
can only occur in the positions g, 7, 1, and o. If 2 occurs in the 0 or | 
position, it can be brought to the j position by one or two applica- 
tions of £. If 2 occurs in the j position, it can be brought to the g posi- 
tion by an application of A. So one can take g=2, m=15 with no loss 
of generality. If 3 occurs in the / or 7 position, it can be brought to the 
0 position by one or two applications of E. Hence one takes o=3, 
e=14, 1=4, c=13. Now d=20—5, and by Theorem 2, /=17—5, 
j=b-—3. By Theorem 1, f=19—b, so that h=b—1, 
p=b—2. Hence we have to find values of b such that b—3, b—2, 
b—1, b, 17—b, 18—b, 19—b, and 20—d are 5, 6, 7, 8, 9, 10, 11, and 12 
in some order. Clearly }=8 or b=12. With )=8, we get 


1 8 | 13 | 12 
14 | 11 2 7 
| 4 5 | 16 9 
15 | 10 3 | 6 


By applying C, B, and A successively to this we get the square that 
would result from putting }=12. Thus all diabolic magic squares can 
be obtained from the one shown by successive applications of A, B, 
C, D, and E. 


THEOREM 4. There are 384 diabolic magic squares of order four. 


Proor. If X and Y are two transformations, denote by X Y the 
transformation effected by first performing Y and then performing X. 
Then (a) BA=AB?*; (b) CA=AD; (c) DA=AC; (d) CB=BD; 
(e) DB=BC?; (f) DC=CD; (g) EA =AB?°CDE; (h) E?A = B°CDE?; 
(i) EB=B*CDE; (j) E*B=ABCDE?; (k) EC=C*E?; (1) E?C 
=B*CD°E; (m) ED=ACE?; (n) E?7D=AB’DE; (0) A?=B*t=C* 
= D4= E*=the identical transformation. These relations were ob- 
tained by inspection and can be easily verified. By means of (a)—(0) 
one can, in any product of A, B, C, D, E, get all the E’s to the right, 
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then all A’s to the left, then all B’s next to the A’s, and finally all 
C’s next to the B’s; so that any product of A, B, C, D, E is equal toa 
product of the form A*B*C7D*E*. Moreover A, B, C, and D all carry 
rows into rows and columns into columns, or rows into columns and 
columns into rows, so that neither E nor E? is a product of A, B, C, 
and D. Moreover, if a square is mapped on a torus, B, C, and D all 
leave the orientation unchanged, so that A is not a product of B, C, 
and D. Clearly B, C, and D are independent. Hence A*B*°C’D*E* 
=A*B’C*D‘E: if and only if a=a (2), B=b (4), y=c (4), 6=d (4), 
and e=e (3). Hence the order of the group generated by A, B, C, D, 
and E is 2X4X4X4X3=384.* As each transformation is a permuta- 
tion of the sixteen elements of the square, different transformations 
cannot yield the same square. So there are 384 diabolic magic squares. 

This theorem disagrees with the result of D. N. Lehmer, A census 
of squares of order 4, magic in rows, columns, and diagonals, this Bulle- 
tin, vol. 39 (1933), p. 981. Lehmer’s count was 48. It is easy to see 
the source of Lehmer’s error. Lehmer found all squares which were 
normalized with respect to certain transformations, R, S, T, and U, 
and picked out the diabolic ones. However, Lehmer’s transformation 
U destroys diabolism, thus certain of his normalized squares which 
are not diabolic give rise to diabolic squares when operated on by U. 
Hence Lehmer missed these squares in his count. 

The present count was checked against Frenicle’s list of magic 
squares. In so doing, we discovered the following startling result: 

If there are two antipodal elements of a magic square of order four 
which add up to 17, then the magic square is diabolic. 

The defining conditions of a magic square plus the single condition 
that a particular pair of antipodal elements add up to 17 give ten in- 
dependent equations in the sixteen elements of the square, as con- 
trasted with twelve independent equations arising from the defining 
conditions for a diabolic magic square. That the two sets of equations 
are equivalent is due to the diophantine condition that the elements 
must be 1, 2, - - - , 16 in some order. 


* H. S. M. Coxeter has suggested that this group may be isomorphic with the 
symmetry group of the Cartesian frame in four dimensions, as the two groups are of 
the same order. That this is so can be shown as follows. The elements Ri=ABCE, 
R.=AB*D, R;=BCE, R,=AB, will generate our group, since A = Ri R3, B= Ri 
E= R,R3R2R;, as may be verified by the relations 
(a)-(0). The R; satisfy the relations R,?=(R,R3)?=(RiR,s)?=(R2Rs)? =(RiR2)* 
= (R2R;)? = (R3R,)4=1, hence our group is a factor group of the above mentioned sym- 
metry group (Coxeter, Journal of the London Mathematical Society, vol. 10 (1935), 
pp. 21-25). The group in question is designated by [32, 4]. Since the two groups have 
the same order they are isomorphic. 
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If one weakens the requirements still further and only asks that a 
square be magic in the rows and columns, then a pair of antipodal 
elements can add up to 17 without the square being diabolic. This is 
illustrated by 


1 15 4 14 
12 2 9 11 
8 7 16 3 
13 10 5 6 


which is magic in rows and columns, but not in diagonals, and which 
has a+k=e+o=17. 

An analogous treatment of the problem of finding all diabolic magic 
squares is given by Kraitchik on page 167 of his book, La Mathé- 
matique des Jeux, where he shows that all diabolic magic squares can 
be derived by successive applications of A, B, C, and D from three 
particular ones which he gives. 


CorRNELL UNIVERSITY 


A NOTE ON REGULAR BANACH SPACES* 
B. J. PETTIS 


Introduction. For an element x of a Banach space Bof it is well 
known that the functional 


X Af) = f(x) 


defined over B,=Bo, the Banach space composed of all linear func- 
tionals (real-valued additive and continuous functions) defined over 
Bo, is linear; moreover{ 


hence the additive operation X,=T7(x) from Bo to Bz =B, is continu- 
ous and norm-preserving. In Bz let B, denote the set of image ele- 


* Presented to the Society, October 30, 1937. 

+ S. Banach, Théorie des Opérations Linéaires, Warsaw, 1932, p. 53. We shall use 
Banach’s terminology. 

¢ Banach, loc. cit., pp. 188-189. 
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ments of T. A space By is regular* if and only if 
(0. 1) Be BL. 


Examples of such spaces are L?, 1”, (1<p< ©), and the hyper-Hilbert 
spaces. 

The purpose of the present note is to list in §1 some conditions 
necessary and sufficient for regularity, to give the form of the general 
linear operation from a regular space to the space /' of absolutely 
convergent series (§2), and to generalize some theorems already 
known concerning the differentiability of functions of bounded varia- 
tion defined to a Banach space (§3). A function Y(R), from the ele- 
mentary figures in a euclidean figure Rp to a B space, is of bounded 
variation, or BV, if and only if }-i-1|| Y(R;)|| is bounded over all finite 
sets of non-overlapping elementary figures Ri, - - -, Rx contained in 
Ro. The |.u.b. of such sums will be denoted by var(Y; Ro). The last 
section includes the result (Theorem 7), analogous to that of J. A. 
Clarkson for uniformly convex spaces and those of Dunford and 
Morse,f that every function BV, from a linear interval to a regular 
B-space, is strongly differentiablef a.e. (almost everywhere).§ 


1. Necessary and sufficient conditions for regularity. We prove the 
following theorem: 


THEOREM 1. Each of the following conditions is necessary and suffi- 
cient that By be regular: 

(1.1) B, is regular; 

(1.2) for linear sets in Bz closure is equivalent to regular closure; 

(1.3) Bo ts weakly complete, and for linear sets in Bz regular closure 
1s equivalent to weak closure as a set of functionals over By; 

(1.4) B. is of the second category in Bo; 

(1.5) each linear functional over By is uniquely extensible over the 
whole of Bz. 


(0.1) implies (1.1). The operation fp = T(F), adjoint (or associate) 


*H. Hahn, Uber lineare Gleichungssysteme in linearen Raéumen, Journal fiir die 
Reine und Angewandte Mathematik, vol. 157 (1927), pp. 214-229. 

t J. A. Clarkson, Uniformly convex spaces, Transactions of this Society, vol. 40 
(1936), pp. 396-414; N. Dunford and A. P. Morse, Remarks on the preceding paper of 
James A. Clarkson, ibid., pp. 415-420. 

t The definition of strong and weak differentiability will be given in §3. 

§ See also Gelfand, Zur Theorie abstrakter Funktionen, Comptes Rendus de 
l’Académie des Sciences de l’URSS, vol. 17 (1937), Theorem 3. This theorem of Gel- 
fand’s was called to the writer’s attention after the present paper had been accepted 
for publication. 
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to T and from B;=B, to By, assigns to each F(X) the functional 
fr(x) defined by 


fe(x) = F(T(«)) = F(X,); 
hence for every Xz in BL 
F(X.) = fr(x) = X.(fr). 


Since B,“ = Be, this implies for every F in B; an fr in B, such that 
F(X) =X (fr) for all X in B.; thus B, is regular. 

(1.1) implies (1.2). If M is a closed linear set in Bz, then for Xo 
not in M there exists an Fo in B; such that Fo(Xo) =1, and Fo(X) =0 
for X in M. Since B, is regular, there is an fo in B, such that Fo(X) 
=X (fo) for all X in Be. Hence X o(fo) =1, X (fo) =0 for X in M, and M 
is a regularly closed set of functionals over B,. 

(1.2) implies (1.3). In general, regular closure implies weak closure 
as a set of functionals, and this, in turn, is stronger than ordinary 
closure. From (1.2) we then have closure, weak closure as a set of 
functionals, and regular closure, all equivalent for linear sets in Be. 
Moreover, the closed linear subspace B.™ is therefore weakly closed 
as a set of functionals over B,; hence By is weakly complete. 

(1.3) implies (0.1). From the weak completeness of Bo it follows 
easily that B, is weakly closed as a set of functionals. From the re- 
mainder of (1.3) the linear set B,. is regularly closed. Hence if there 
existed an Xo in B, that was not in B,, there would then be an fo 
such that Xo(fo) =1, while X(f.) =0 for all X in BX. But then we 
would have fo(x) =0 for all x in Bo, that is, fo would be the identically 
zero functional over Bo, contrary to Xo(fo) =1. 

That (1.4) is equivalent to (0.1) follows from the fact that the set 
of image points of a linear operation from one Banach space to an- 
other is either of the first category in the contradomain space of the 
operation or consists of the whole of that space. 

If B. contains an Xo not in B,., then, since B. is linear and 
closed, there is an Fy in B; such that Fo(Xo) =1, and Fo(X) =0 for 
X in B, . Thus the zero functional over B. is extensible in at least 
two different ways, which contradicts (1.5). 


THEOREM 2. If Bois separable, then each of the following conditions is 
necessary and sufficient that Bo be regular: 
(1.6)* Bois weakly complete and weakly compact; 


* (1.6) is well known; for the sufficiency see Banach, loc. cit., p. 189, and for the 
necessity, T. H. Hildebrandt, Linear functional transformations in general spaces, this 
Bulletin, vol. 37 (1931), p. 200. 
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(1.7) Bo is weakly complete and B, is separable; 
(1.8) im By, weak convergence as a set of functionals is equivalent to 
weak convergence to an element as a set of elements. 


If Bo is separable, so is Bx ; hence if Bo is in addition regular, then 
B:, and therefore B,, are separable. Thus from (1.3) the regularity 
and separability of By together imply (1.7). 

Condition (1.6) follows from (1.7). If {xn} is a bounded sequence 
in Bo, then, by diagonalizing, a subsequence {x,’ } can be found such 
that lim,f;(x,’) exists for every f; in a denumerable set dense in B;. 
The bounded sequence {x,’ } is then weakly convergent in Bo. 

Thus if Bo is separable, conditions (0.1), (1.6), and (1.7) are all 
equivalent. 

The necessity of (1.8) is clear. If lim,f,(x) exists for every x, then 
fo(x) =lim,f,(x) is a linear functional over Bo, and if Bo is regular, 
then F(f,,) converges to F(fo) for every F in Be. If (1.8) holds and F(f) 
is linear over B,, then the convergence of {f,} to fo, as a sequence of 
functionals over Bo, implies lim F(f,) = F(fo); hence F(f) is weakly 
continuous over B;. Since By is separable, this insures the inclusion of 
F in BL 


THEOREM 3. Jf Bo is regular, so is every closed linear subspace By . 


Suppose X’(f’) is a linear functional over the space By conjugate 
to By. Define X(f) = X’(f’), where f varies over Bi, and f’ is the linear 
functional over By’ determined by f. Since Bo is regular, there exists 
an x in By such that X(f) =f(x) for every f in B;. We have only to 
show that x is in By’; for then 

X"(f’) = X(f) = f(*) = f'(*) 
for all f and hence for all f’, since any f’ can be extended to form an f. 
If x is not in the closed linear subset By’ of Bo, there exists an f* in B, 
such that f*(x) ~0, f*(x’) =0 for all x’ in By’. But then X(f*) =f*(x) 
~0, while X (f*) = X’(f*’) = X’(0) =0, a contradiction that completes 
the proof. 

If By is regular, then any separable closed linear subspace is also 
regular and is therefore weakly complete and weakly compact. The 
space By must then have these two properties also. This gives us the 
following corollary: 


COROLLARY. A regular space is weakly complete and weakly compact. 


2. Operations from a regular space to the space /'. We prove the 
following theorem: 


{ Banach, loc. cit., p. 131. 
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THEOREM 4. Jf Bois regular and the operation 
T(f) = {f'} 
from B, to l' (space of absolutely convergent series) is linear, then T is 
completely continuous and is defined by a series >. *x,, unconditionally 


convergentt in Bo; conversely, such a series defines a linear completely 
continuous operation from B, tol’. The norm of 


Given a linear 7, let F= T(m) be the adjoint linear operation from 
l'=[®, the space of bounded sequences, to Bs. If F; = 7(m;), where m; 
is the linear functional over /' defined by the sequence having 1 in 
the ith place and 0’s elsewhere, then 


Fi(f) = m(T(f)) = mA {f*}) = 
Since Bo is regular, there exists an x; such that 
fi = PAD = f(x) 
for all f in B,. Hence 
(2.1) T(f) = {f(x}. 


From the weak completeness of By and the fact that Deel f(xs)| <0 
for all f, the series }>x; is unconditionally convergent;{ moreover, 
given e>0O there is an N, such that§ 


(2.2) | Mx) | <e 

Ne+1 
for all f with ||f|] <2. If ||f,*|| <1, then by diagonalizing we can find a 
subsequence {f;} such that lim:f,(x,) exists for every n. This fact, 
combined with (2.2), gives 


| — f | < 2e 


n=1 


for i, 7 greater than a suitably chosen K,. Thus 7(f) is completely 
continuous. 

Conversely, suppose >.x, converges unconditionally in By. The op- 
eration 7(f) in (2.1) is then defined and additive from B, to /'. Since 


+ A series is unconditionally convergent if every reordering is convergent; 
W. Orlicz, Beitrége zur Theorie der Orthogonalentwicklungen, 11, Studia Mathematica, 
vol. 1 (1929), pp. 241-255; p. 242. For a statement of equivalent definitions see 
Banach, loc. cit., p. 240. 
t Orlicz, loc. cit., Theorem 2. 
§ Orlicz, loc. cit., p. 246. 
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it is the “pointwise” limit of the linear operations 
(2.3) T,(f) = {f(x1),--- , fxn), 0, 0, -- n= 1, 2,°°° 


T must also be linear. The complete continuity follows from the first 
part of the proof. 

The statement concerning the norm of T is an immediate conse- 
quence of (2.1). 

Without resorting to (2.2) the complete continuity of T could have 
been proved by using the following facts: (i) T is linear, (ii) B; is 
weakly complete and weakly compact, (iii) the property of converg- 
ing weakly to an element is preserved under linear operations, (iv) in 
l' weak and strong convergence are equivalent. 

In view of (1.1), Theorem 4 can be stated as follows: 


THEOREM 4’. Jf Bo is regular, an operation 
T(x) = 


from By to l' is linear if and only if it is defined by an unconditionally 
convergent series > f;, in By. If it is linear, it is completely continuous and 


| = Lub. |. 


It is clear from the forms of their general linear functionals that the 
spaces L?, 1?, (1<p< ©), are regular. Applying Theorem 4’ and (2.2) 
to these spaces we obtain the following corollary: 


COROLLARY. f An operation T(x) from L? (l*), (1<p<~), tol‘ is 
linear if and only if there exists in L”’ (l’), (p’=p/p—1), @ sequence 
{ys} = = {fat} }) such that 


(1) >| f yi(s)x(s)ds | < 


i=1 
for every x=x(s) (x= {xi}) in L? that is, converges uncondi- 
tionally in L”’ (1”’), and 


m= 4f (7) = { weit), 


If T is linear, then 


j=1 


{ In this corollary the case /? to /' is already known. See H. R. Pitt, A note on 
bilinear forms, Journal of the London Mathematical Society, vol. 11 (1936), pp. 174- 
180; and L. W. Cohen and N. Dunford, Transformations on sequence spaces, Duke 
Mathematical Journal, vol. 3 (1937), pp. 689-701. 
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(i) T is completely continuous; 
(ii) given €>O0 there exists an N, such that 


f 


x 


i=Ne¢ 


> 


j=l 


uniformly for all x with 


fi x(s) |"ds (2 wal” < ); 


j=1 


f 


j=1 | 


x 


- l.u.b. > 


i=1 


(iii) 7 


3. Differentiation in regular spaces. In this section we consider the 
differentiability of functions BV from a linear interval to a regular 
space. Our main purpose is to prove the theorem stated in the intro- 
duction. Every such function is strongly differentiable a.e. A function 
Y(R), from the elementary figures in a euclidean figure Ro to a 
Banach space Bo, is said to be strongly differentiable at a point s in 
the interior of Ro if there exists an element y(s) of Bo such that for 
cubes J lying in Ry and containing s it is true that | Y(1)/| 1| —y(s)|| 
tends to 0 with |Z, the measure of J. The function Y(R) is weakly 
differentiable at s if there exists a y(s) such that 

lim | — | = 0 


sel 


for every linear functional f. 


THEOREM 5. If Y(R) is defined and BV from the figures in a euclidean 
figure Ry toa Banach space Bo, and if Y(R) is weakly differentiable a.e. 
in Ro to the function y(s), then y(s) is Bochner integrable.* 


Let J, be a cube containing Ro, and {II,,} a sequence of partitions 
of J, into a finite number of non-overlapping (except on their bound- 
aries) and non-degenerate subcubes, with II,; a repartition of II, and 
lim,.. norm II, =0. If s lies in the interior of a cube /;. of II, and Ro 
contains J;,, then define y,(s) as Y(Iin)/| Tin| ; otherwise let y,(s) van- 
ish. Then a.e. in Ry we have {y,(s)} converging weakly to y(s), hence 


*S. Bochner, Integration von Funktionen, deren Werte die Elemente eines Vektor- 
rdéumes sind, Fundamenta Mathematicae, vol. 20 (1933), pp. 262-276. 
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lim inf || 
a.e. in Ro. From Fatou’s lemma and the inequality 
f || ya(s)||ds S var (Y; Ro) < ~, 
Ro 


the real-valued lim inf,||y,(s)|] must be summable. The abstract func- 
tion y(s) is measurable in the sense of Bochner,* since a.e. in Ro it is a 
weak derivative. Hence ||y(s)|] is measurable and summable, which 
implies that the measurable y(s) is Bochner integrable. 


Coro.iary.{ If, in addition to the assumptions in Theorem 5, Y(R) 
ts additive and AC (absolutely continuous), then 


(A) = J soa, 


and Y(R) is strongly differentiable a.e. to the function y(s). 


In view of the strong differentiability a.e. of the Bochner integral,§ 
it is sufficient to prove (A). Let Z(R)=/J,y(s)ds; then|| f(Z(R)) 
= for every f in = By and every Rin Ro. But f(Y(R)) is 
a real-valued additive and AC function that is differentiable a.e. to 
f(y(s)), and hence f(¥(R)) for 
every f and R. Thus Y(R) =Z(R) for every R. 

Because of this corollary a theorem of Dunford and Morse] may 
be altered to read: 


THEOREM 6. A Banach space By has the property (DBV), namely, 
that every BV function from a linear interval to By is strongly differenti- 
able a.e., and its derivative 1s Bochner integrable, if and only if every 
Lipschitzian function from a linear interval to By is weakly differenti- 
able a.e. 


THEOREM 7. Jf By is regular, then Bo has property (DBV). 


* Bochner, loc. cit., p. 263. 

t See On integration in vector spaces, a paper by the present author which is to ap- 
pear in the Transactions of this Society. 

t Theorem 5 and the corollary are generalizations of theorems due to Clarkson, 
loc. cit., pp. 409-410. 

§ Bochner, loc. cit., p. 269. 

|| N. Dunford, Integration and linear operations, Transactions of this Society, vol. 
40 (1936), pp. 474-494; p. 475, Theorem 2.3. 

Loc. cit., p. 415. 
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Suppose Y(¢) is defined and Lipschitzian from a linear interval 
[a, 6] to a regular space Bo. From the Lipschitzian property it is 
clear that every value of Y(t) lies in the separable closed linear sub- 
space By generated by the values of Y(¢) at the rationals. By (1.7) 
and Theorem 3 there exists a denumerable set {fi } dense in the space 
By adjoint to Bo ; and By’ is weakly complete. Let Y(R) be the addi- 
tive figure function determined by Y(t). Then f/(Y(R)) is a real- 
valued additive Lipschitzian figure function; hence there exists in 
[a, b] a measurable set S with |S| =b—a and such that s in S im- 
plies f/ (Y(R)) differentiable at s for each 7. Thus if s is in S, and if 
{Z,} are intervals closing down on s, it follows that 


exists for every f’ ina set dense in B, ; moveover, by the Lipschitzian 
property | Y(I,)/ | Zn] | | <K. These two conditions imply that 
lim f’( YC, )/ | In|) exists for every f’ in Bj, that is, { VU, )/|In | } 
is a weakly convergent sequence. Since B,’ is weakly complete there 
is an element y(s) to which this sequence converges weakly. If {I, } is 
another sequence of intervals closing down on s, then { (Yd, )/| ER | } 
is also weakly convergent, and for every f’ 


lim f’(¥(In)/| In|) = lim In|) = 


Thus Y(R), as a function defined to By, is weakly differentiable 
a.e. in [a, bj. Since any f in B, defines an f’ in B/ , Y(R), as defined to 
Bo, has the same property. On application of Theorem 6, the proof is 
complete. 

If L denotes the space of functions Lebesgue integrable on [0, 1], 
we have the following corollary: 

Coro_iary. If By is regular, then T(@) from L to Bo is linear tf and 
only if there exists an essentially bounded and Bochner measurable func- 
tion x(s) on [0, 1 ] to By such that 


(6) = J 


the integration being in the Bochner sense. The norm of T ts ess., sup. 
x(s)|| 
To obtain this from the theorem the reader is referred to a proof 
by Dunford.* 
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RELATIONS AMONG THE FUNDAMENTAL SOLUTIONS 
OF THE GENERALIZED HYPERGEOMETRIC 
EQUATION WHEN /=¢+1. 

I. NON-LOGARITHMIC CASES* 


F. C. SMITH 


It is the purpose of this paper to develop the relations among the 
non-logarithmic solutions of the equation 


Pp q 
(1) {Tl @ +0) - @+a- by =o, p=qt1, 
t=1 
where @ =2(d/dz). The results of this paper generalize those of Mehlen- 
bacher,t who studied the case in which p=2 and g=1, and extend 
those of Barnest who obtained the asymptotic developments of the 
non-logarithmic solutions of (1) for those cases in which p<q+1. 
The succeeding analysis will be simplified if we rewrite equation 
(1) in the eq uivalent form 


q+1 1 q+1 
(2) Con = 1. 
t=1 t=1 


If no two of the a; or c; are equal or differ by an integer, then equation 
(2) has g+1 linearly independent solutions about the point z=0, 
which may be written 


Td +o — + a — + 2) 
1; <1; 


and g+1 linearly independent solutions about the point z= © which 
may be written 


(4) wnt — + Gj) T(1 — + 03+) 2* 

j =1,2,3,---,q +13 cer = 1; > 1. 


* Presented to the Society, December 30, 1937. The logarithmic cases will be con- 
sidered in a later paper. 

+ See this Bulletin, Abstract 42-5-169. 

t E. W. Barnes, The asymptotic expansions of integral functions defined by general- 
ized hypergeometric series, Proceedings of the London Mathematical Society, (2), vol. 
5 (1907), pp. 59-116. 
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In this paper, we shall develop the relations among solutions (3) and 
(4) of equation (2). 
In this connection, we may state the following theorem: 


THEOREM 1. If no two of the a, or c, are equal or differ by an integer, 
then the solutions (3) of equation (2), when extended analytically outside 
their circle of convergence, may be expressed linearly in terms of solu- 
tions (4) in the following forms, it being understood throughout that 
0 <argz<27: 


Yo; = — ax)P(1 + ce — 63) 
k=1 ax) 
(5) il T(a, — a, + — ¢;) 


j= 1,2,3,---,¢t¢ = 1.* 


ProoF. It is well known that the function 


+ a, — + w) 
6 y= 
(6) g(w) 


satisfies the relation 


g(w) = [1 + y(w)], lim y(w) = 0, | arg w | 


Therefore, 


| e(w) | < | 


(at—ct) 


| [1 + | v(w) | ] = + | |], 


where R and J denote the real and imaginary parts of 5 as—c,), 
respectively, and where ¢ =arg w. For sufficiently large values of | w| : 
we have | y(w) | <e, larg w| <m. Thus, for such values of | 2 | 


| g(w) | < Re-I¢(1 + €) < | w *(1 K|w|, 


where K is a constant independent of w. The function g(w), therefore, 
satisfies the conditions of a general theorem due to Ford,{ but, as ap- 


* As far as relation (5) itself is concerned, the restriction cg,:=1 is unnecessary. 
The restriction is necessary, however, in order that the Yo; and the Y,,; be solutions 
of equation (2). 

+ See, for example, Milne-Thomson, The Calculus of Finite Differences, Macmillan, 
1933, p. 255. 

t See Ford, The Asymptotic Developments of Functions Defined by Maclaurin Se- 
ries, University of Michigan Press, 1936, pp. 5, 9. 
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pears from (6), has simple poles at the points 
(7) w=c;—a,—n—1, n=0,1,2,---;k =1,2,3,---,qg4+1. 
An application of this theorem to the function 


(8) {@ = = XT] (1+ a,—¢;+n) 


n=0 n=0 t=1 Cj + n) 


gives the result 


1 ert 


(9) ami 2" k=l n=0 
0 < arg z < 27, 


where R,,,, represents the residue of the function 


a(— 2)"g(w) + a,—c;+ w) 


at the point w=c;—a,—n—1. But the first series on the right in (9) 
drops out since the denominator of every term contains a factor for 
which ¢t=j, and 
1 1 
T(i+ec;—c;— 72) — n) x 


0, 


It follows that our problem reduces to that of computing the residue 
term. 
For this purpose, consider 


1 + a, + w) 
-—f — z)*T(w)T(1 — di 
(1) Ra = 55 J, 


where C,,,, surrounds the pole w=c;—a,—n—1 of (10) but no other 
pole of (10). If in (11) we replace w by —(w—c;+a,+n+1), we ob- 
tain 


( n—1 


2ri 
2r(—w— 


s+! T(a,—a,—w—n) 


dw, 
w—n) 
t¥k 
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where Cy surrounds the origin. By several applications of the relation 
T 

(13) r(i — w) = 

we may change (12) to the form 


(- — 1)* 


Raye = — 
2xi 
— — — + a + 1) 
sinww — a, — — + + 1) 
(14) T(a, — a, — w)T(w — + ax + 1) 
— ay — w)T(w — ce + ax + 1) 


Mw — +1) 
T(w— a +a, +n-+ i) 


which by virtue of a well known theorem in the calculus of residues 
is equal to 


* 


(— — (1 — + ax) 


= 
(15) 2" — ax)T (1 — ce + ax) 
T(a, — ax)T (A — ar) +1 m) 
wma — ax) (1 — + m1 n) 
Therefore, 
— T'(c; — (1 — c; + ax) 
n=0 — (1 — + ax) 


+1 -—-atatn)i 


(16) 


T(c; — ax)T(1 — + ax) 


(cx ax) 
q+1 — ar) = ax) 
- Y 
t=1 T'(c, — ox) ax) 
t#k 


" See Ford, op. cit., p. 5, equation (20). 


Thus, 
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(17) pers k=l — ax) 
qt+1 T(a, a) 
tk 


Substituting (17) into (9) and recalling that the first series of (9) is no 
longer present, we have 


+1 
f@) ~ — ax)T(1 — + ar) 
ax) 


(18) 
tm1 — ax) 
Since each Y,,, which appears in (18) is defined by a series which con- 
verges for | z| >1, we may replace ~ by = in (18). Then, upon noting 
that 
ott + — 
(19) Yo; = 
we see that (18) is equivalent to the desired result (5). 
By means of a similar proof, we may establish the following theo- 
rem: 


THEOREM 2. If no two of the a; or c; are equal or differ by an integer, 
then the solutions (4) of equation (2), when extended analytically outside 
their circle of convergence, may be expressed linearly in terms of solu- 
tions (3) in the following forms, it being understood throughout that 
0<arg (1/z) <2r7: 


+1 
= a;)T(1 a. + a;) 
= ax) 


k=1 
(20) — — ae + aj) 
(cx — — + aj) 


t#k 


1, 2,.3,°+: Ces = 1. 
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A NOTE ON FREDHOLM-STIELTJES 
INTEGRAL EQUATIONS* 


F. G. DRESSEL 


1. Introduction. The object of this paper is to show that the inte- 
gral equation 


(1) fla) = m(x) +r f 9), 


IA 


1, 


can be changed into an ordinary Fredholm equation when G(x, y) is 
absolutely continuous g(y).{ The integration is carried out in the 
Young-Stieltjes sense, and g(y) is a bounded, monotone increasing 
function. 


2. Lemmas. If h(x) is of bounded variation and we set h(x) =h(0), 
(x<0), and h(x) =h(1), (x>1), then we may define the completely 
additive function of sets h(e) by 

h(e) = h(x. + 0) — h(x, — 0), e= e(%, St S x). 


Using this notation we have the following lemma: 


Lemma 1. If f(x) is measurable Borel then 


the left side being Young-Stieltjes integration, the right Radon-Stieltjes. 


In case one integral does not exist the equality sign is taken to 
mean that the other integration is non-existent. Because of the prop- 
erties of the integrals under consideration, we need only prove the 
equality for the functions 


f(t)=1,2=0, fils) 
= 0, 


0, elsewhere. 


* Presented to the Society, December 29, 1936. 

+ For a discussion of (1) see G. C. Evans and O. Veblen, The Cambridge Colloquium 
Lectures on Mathematics, American Mathematical Society Colloquium Publications, 
vol. 5, 1922, p. 101. 

t For terminology see Alfred J. Maria, Generalized derivatives, Annals of Mathe- 
matics, vol. 28 (1926-1927), pp. 419-432. I am much indebted to Mr. Maria for 
many valuable suggestions. 

All functions used in the present paper are assumed to be measurable Borel. 
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We have 


f fi(x)dh(x) = h(a + 0) — h(a — 0) = h(a) -f fi(x)dh; 


f fo(x)dh(x) = h(B — 0) — h(a + 0) = Ae) -f fo(x)dh, 


where is the open set a<t<£.* 


Lemna 2. If G(x) is absolutely continuous with respect to the bounded 
monotone increasing function g(x), then 


f = f 


where DG(x) is the derivative or one of the derived numbers of G(x) with 
respect to g(x). 


Mr. Mariat has made the important step in the proof of the lemma 
by showing that 


+ 0) — — 0) = f DG(x)dZ, 
E 


where E is the set x; ¢<-x2. For the function f,(x), making use of 
Lemma 1, we have 


f = Ge + 0) Ge 0), 


f fi(x)DG(x)dg(x) -{ fi(x)DG(x)dz -f DG(x)dz 


= G(a + 0) — G(a — 0), 


where E is the point a. For f2(x) we have, if e is the open set a<x<B, 


f folx)dG(x) = GB — 0) — G(a +0), 


= G(6 — 0) —G(a + 0). 


From the above material the lemma readily follows. 


* The same reasoning shows that /¢f(t)dh(t) is equal to /of(t)dh, for 0<xX\, if 
h(t) is continuous from the right except perhaps at x=0. 
Tt Loc. cit., p. 430. 
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3. Transformations. Our first theorem is the following. 


THEOREM 1. Jf G(x, y) ts absolutely continuous g(y) then equation (1) 
can be written in the form 


(2) f(x) = m(x) +2 f K(x, »)f(y)dg(y), 


where K(x, y) =DG(x, y), the derivative being taken with respect to g(y), 
a bounded monotone increasing function. 


This is immediate from Lemma 2. 


THEOREM 2. If m(x) and K(x, y) are bounded, then the solution of (1) 
and (2), except for characteristic values of d, can be written 
f(x) = m(x Dn) m(y)dg(y), 


where 


D(A) = 1 s)dg(s) +---, 
0 


1 


K(x, y) K(x, s) 
K(s, y) K(s, s) 
The proof follows along the same lines as in the ordinary case. We 


now state a corollary of Theorem 2 that represents most of the known 
results concerning solutions of equation (1). 


Coro.iary.* | G(x, y2) —G(x, | —g(y)|, then, except- 
ing characteristic values, equation (1) has (3) as a solution. 


dg(s) +---. 


D(x, ¥;) = K(x, 9) — f 


0 


Any result for the ordinary Fredholm equation carries a related re- 
sult for equation (1). To see this, we assume without loss of generality 
that g(y:) <g(y2) if and apply to (2) the transformation ft 


B(s) = lim sup E,(s 2 g(y)), g(0) Ss S g(t), 


(4) f(z) = m(x) + K(x, y)f(y)dg(y) 


g(1) 
m(x) K(x, 
9(0) 

* This includes the case handled by W. C. Randels, On Volterra-Stieltjes integral 
equations, Duke Mathematical Journal, vol. 1 (1935), pp. 538-542. 

¢ Banach, Théorie des Opérations Linéaires, Warsaw, 1932, p. 6. 


= 
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If we let w be any of the possible solutions of 
= Bw), g(0) gil), 


we may write (4) in the form 
g(1) 


F(w) = M(w) + af k(w, s)F(s)ds, 


g(0) 


where Fw) =f(8@)), M@)=m(B@)), k@, s)=K(B(»), B(s)). We 


thus have our main result: 


THEOREM 3. When G(x, y) is absolutely continuous g(y) the Fred- 
holm-Stieltjes integral equation (1) is reducible to an ordinary Fredholm 
integral equation. 


4. Mixed linear equations. The mixed equation* 
m 1 
(5) f(x) = m(x) + (x) +d f K(x, s)f(s)ds 
i=l 0 


can easily be put into the form 
1 
f(x) = m(x) + rf R(x, s)f(s)dg(s). 
0 


Thus from Theorem 3 we see that equation (5) is reducible to a 
Fredholm integral equation. 


DuKE UNIVERSITY 


A THEOREM ON QUADRATIC FORMS?{ 
WILLIAM T. REID 
In this note the following result is proved: 


THEOREM. Suppose A [x]=a.sxaxs,t Bl[x]=Dasxaxg are real quad- 
ratic forms in (xa), (a=1,---+, m), and that A[x]>0 for all real 
(xa) #(0a) satisfying B[x]=0. When there exists a real constant Xo such 
that A [x]—oB[x] is a positive definite quadratic form. 


This theorem is of use in considering the Clebsch condition for 
multiple integrals in the calculus o/ variations. A. A. Albert§ has given 


* W. A. Hurwitz, Mixed linear integral equations of the first order, Transactions 
of this Society, vol. 16 (1915), pp. 121-133. 

Tt Presented to the Society, December 30, 1937. 

t The tensor analysis summation convention is used throughout. 

§ This Bulletin, vol. 44 (1938), pp. 250-253. 
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an algebraic proof of this result. The following proof is more analytic 
in character than that of Albert and is of interest because in the proof 
itself we obtain directly, in terms of the roots of the characteristic 
equation (1) below, the interval in which the value A, of the theorem 
may be chosen. 

Suppose that the matrix ||b,s|| has rank n—r, (0Sr<n), and let 
Ua =Uaox, (K=1, - - - ,7), be linearly independent solutions of the equa- 
tions bagug=0, (a=1,---, 2). 


Lema 1. The characteristic equation 
(1) D(A) = | aas — Abas| = 0 


is a polynomial of degree n—r in X. 


Let (S=1,---, be sets linearly independent of 
the above defined u.. and such that the determinant | was|, 
(a, B=1,---, 2), is different from zero. Clearly the degree of 


D(A) is equal to that of the polynomial D,(A) =| (dey—Nbay) tty] « 
Since the first r columns of D,(A) are independent of A, we obvi- 
ously have d*D,/dX\*=0 if k>n-—r. Also, the (n—r)th derivative is 
equal to (n—1)!| — 
and if its value is zero, there are constants (¢.)#(0.) such that 
Ue Ua =UarrsCris Satisfy the equations =bagug. By use 
of these equations together with u2b.3=03 we have A [u°]=0; hence 
(u£)=(0.), by the hypothesis of the theorem, and consequently 
(c,) = (0,). This implies in turn that =0.; and since ) 
and is independent of the sets wa., (k=1, - - - , 7), we have a contra- 
diction to the assumption that || bag|| is of rank —r. The above lemma 
is therefore proved. 

Corresponding to a root of (1) the number of linearly independent 
solutions of the associated equations 


(2) (das Abas) = a, B = Nn, 


is termed the index of \ asa root of (1). For sets x.’ =(x2' ), xa’ =(x2'") 
we denote by A[x’; x’’], B[x’; x’’] the bilinear forms xs’, 
DapXa x3’. The first of the following lemmas is immediate. 


Lemma 2. If \’, X”’ are distinct zeros of D(X), and ya , ya’ are solu- 
tions of the corresponding equations (2), then B|y’; y"’]=0=A [y’; y’’]. 


Lema 3. The zeros of D(d) are all real, and the corresponding solu- 
tions of (2) may be chosen real. 
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For suppose \’ =A, +72, (A2 #0), is a zero of D(A), and ya = Yai 
is a corresponding solution of (2). Then the conjugate values 7, sat- 
isfy (2) for X=)’. By Lemma 2 we have 0=4A [y; §] =A [y:]+A [ye], 
0=B[y; Now if either B[y,] or Bly2] is zero, so 
also is the other, and, in view of the hypothesis of the theorem, 
we have a contradiction. On the other hand, if B[y2]0, let c be 
a value opposite in sign to \2/B[y,] and such that 0=B[y,+cye] 
= B[y:]+2cB[y1; y2]+c?B [ye] =(1—c*)B[y: ]+2cB[y1; Using the 
relations satisfied by ya, Yaz One obtains, by direct calculation, 
A [yi +cy2] = (2che/Bly])({ Bl] y2]}?)<0, which is a 
contradiction. Hence there are no complex zeros of D(A), and the 
solutions of equations (2) may be chosen real. 


LemMA 4. If X\=X’ is a zero of D(A), then its index is equal to its 
multiplicity. 


Suppose the, index of X’ is equal to k, and denote by ya=Yaj, 


(j=1,---, k), corresponding linearly independent solutions of (2). 
Now choose yan, (2=k+1, - - - , m), such that the determinant | ; 
(a, B=1,---,), is different from zero. Clearly the multiplicity of 


\’ as a zero of D(A) is equal to its multiplicity as zero of 


AA) =| (Gay —Abar) - 
Since the first k columns of A(A) are zero for \=X’, the multiplicity 
of \’ as a zero of A(A) is not less than k. Moreover, 


d*¥A/dd* = k!| — base; (Gas — Yon |, 


3 
(3) j=l,---, =k+1,---,m. 


Now if this derivative is zero, there are constants (c,.) ~(0.) such that 
the set y2 =Yajlj, Va =YanCn Satisfies 


Since (y2) satisfies (2) for \=)’ we find from (4) that B[y®]=0, and 
consequently, using again equations (2), that A [y®]=0. In view of 
the hypothesis of the theorem it then follows that (v2) =(0.) and 
(c;) =(0;). But from equations (4) we see that (y,.’ ) is then a solution 
of (2) for \=X’; and, since this set is independent of the & sets ya;, 
we have a contradiction to the assumption that & is the index of 0’. 
Hence the expression (3) is not zero, and the multiplicity of \’ as a 
zero of D(A) is also k. 

Now let\i, SA2S - - - SAn_, denote the zeros of D(A) each repeated 
a number of times equal to its index, and let ya =Yas, (S=1,---,"—?r), 
be corresponding linearly independent real solutions of equations (2). 
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Lemna 5. The solutions ya, are such that B[y,|#0, (s=1, - --,n—r); 
moreover, these solutions may be so chosen that B[y,; y,| =0 if s¥t. 


The first part of the lemma follows readily from the hypothesis of 
the theorem, since A [y,]=\,B[y.], (s=1, - - -, If the 
second part of the lemma is true by Lemma 2. If 4, =;, the corre- 
sponding solutions may be so orthogonalized in view of the relation 
Bly, | #0. 


Coro.iary. The determinant | uaxYas| is different from zero. 


Lemna 6. If for a value s we have B[y,]>0, then for each t such that 
A, ZA, we also have B[y,|>0. 


For suppose \,=A, and B[y,] <0. If c?=—B[y,]/B[y.], then by 
the above lemma we have (ya:+Cyas)~ (02). Moreover, B[y.+cy, | 
=B[y.]+cB[y.]=0, A [y+ey]=A [y.] 0, 
which is impossible in view of the hypothesis of the theorem. 

Now let A denote the interval of the \ axis defined by \<A, if 
B[y,]>0,A>A,_- if Bl yn_-] <0, if Bly,] <0, >0. 
The interval A is uniquely determined in view of Lemma 6. We shall 
now show that A[x]—X .B[x] is a positive definite quadratic form 
for X>) on A. If (xa) (02), it follows from the corollary to Lemma 5 
that there are unique values (c,, d,)~(0., 0,) such that 
+Yad,. Moreover, since tecbas=0, we have A[u,; y.]=0, («=1, 

-,7;s=1,---,n—r), and A[u,, u,|c.c,20, with the inequality 
sign holding unless (c,) = (0,). Now if Xo is on the above defined inter- 
val A, each of the terms (A,—\o)B[y.] is positive. From the relations 


A[x] — roB[x] = A [ue: + (Aly; — 


= A[uc; + (As — Ao) Bly. ]d? 
s=1 
we therefore conclude that A [x ]—X,B[x] is positive definite for \ on 
A. In conclusion it is to be noted that if \» is not on this interval, the 
form A[x]—AoB[x] is negative or zero for at least one of the sets 
(xa) = (Yas); (s=1, n—r). 


UNIVERSITY OF CHICAGO 
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ON THE EXTENDED FORM OF CAUCHY’S CON- 
DENSATION TEST FOR THE CONVER- 
GENCE OF INFINITE SERIES 


R. E. MORITZ 


Cauchy’s celebrated theorem,* which forms the basis of some of 
the most powerful chain tests for the convergence and divergence of 
infinite series, notably the logarithmic scale and DeMorgan and 
Bertrand’s T and 7 scales, is as follows: 


Let f(n) represent any positive continuous and single-valued function 
of n which is monotone decreasing as n increases. Suppose, moreover, 
that lim, ..f(") =0. Then >,°f(n) and are both convergent or 
both divergent, provided a is a positive integer not less than 2. 


The restriction a <2 was removed by Kohnt in 1882. The proof 
covers several pages. Kohn showed that the theorem holds for any 
real value of a greater than unity. 

A shorter proof of Kohn’s extension of Cauchy’s theorem was given 
by Hillf in 1896. Hill’s proof depends on two auxiliary theorems 
which, in turn, are based on the two previously established inequalities 


1 1 
(1 + + —) < test t,, 
a 


1 
— & 
a’ 


where ¢, is the integral part of a’. 

In view of the fundamental importance of Cauchy’s theorem the 
following simple analytical proof of the theorem with Kohn’s exten- 
sion may, therefore, be of some interest. The proof which follows is 
in fact in the nature of a corollary to the well known integral test 
which may be stated as follows: 


Under the assumptions set forth in Cauchy's theorem, >.f(n) and 
J°f(x)dx converge and diverge together. 


Set x =a", (a>1); then log a>0 and 


* Cauchy, Anciens Exercises, vol. 2 (1827), p. 221. 
+ Kohn, G., Grunert Archiv, vol. 67 (1882), pp. 63-95. 
t Hill, M. J. M., Messenger of Mathematics, (2), vol. 26 (1896-1897), p. 102. 
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sea = log af = log af 


Therefore >. °f(n) converges and diverges with [“af(a*)dx, and it re- 
mains only to show that [“a*f(a*)dx converges and diverges with 

For values of x between the two integers m and +1, we have, if a 
is greater than unity, a*<a*<a"*'; and since f(x) is monotone de- 
creasing, 

f(a") < f@) < fe"). 
Furthermore, since f(x) is positive we have, on multiplying these in- 
equalities term by term, 


a*f(a*t") < a*f(a*) < a"*¥f(a"), 
from which, on integrating with respect to x between the limits n 
and +1, we have 
n+1 


1 
< a*f(a*)dx < a-a"f(a") 
a 


n 


and therefore 

n+2 n+3 

< f a*f(a Me + a*f(a*)dx 


These two inequalities make it obvious that) ,“a*f(a") and [“a*f(a*)dx 
converge and diverge together. 

This completes the proof of both Cauchy’s theorem and Kohn’s 
extension. 
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ON THE BOUNDARY CONDITION 0u/dn+au=0 
FOR HARMONIC FUNCTIONS* 


HILLEL PORITSKY 


1. Introduction. In a recent paper (referred to below as I)f the au- 
thor considered the nature of the “reflection” or analytic continuation 
of harmonic functions across a plane over which the boundary condi- 
tion 


Ou 
(1) —+au=0, a = const., 

On 
applies. In the following this investigation is continued. First, the 
case of a spherical (circular in two dimensions) boundary, over which 
the boundary condition (1) holds is considered. It is shown that sin- 
gularities admit of a similar “reflection,” as in the case of a plane in I; 
thus a point singularity at Py outside a sphere S over which (1) holds 
is reflected into a point singularity at P, (the spherical inverse of Po 
in S) and into a distribution of singularities along the radius vector 
from the center 0 to P. 

Returning to plane boundaries, we apply boundary conditions of 
the form (1) over two parallel planes. A rather complicated “reflec- 
tion” of singularities results, consisting of point singularities as well 
as of distributed line singularities. The point singularities are located 
at the periodic row of points obtained by reflecting the original singu- 
larity Po first in one plane, then in the other one; reflecting these 
images in the two planes; and so forth. The line singularities are dis- 
tributed over the straight line through the above point singularities. 
The density of the distributed line singularities is an analytic function 
of the distance along the line bearing them between the point singu- 
larities, but changes abruptly from one analytic function to another 
one in crossing these points. 

Some of the above features are believed to be typical of analytic 
continuations of a great variety of expansions in characteristic func- 
tions related to two point problems. 


2. Circular and spherical boundaries. We shall consider the ques- 
tion of reflections of singularities of harmonic functions across spheri- 
cal or circular boundaries corresponding to the condition (1). 


* Presented to the Society, September 5, 1934. 
{ This Bulletin, vol. 43 (1937), pp. 873-885. 
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First we consider the two-dimensional harmonic case. Suppose that 
the function x satisfies 
Ou 
(2) —-—au=0, for r=1, 
or 
and is harmonic outside of this circle except near a point Po where it 
becomes infinite like —In7ro; u is also to vanish at infinity. The point 
P, is along the x axis at x=h =e” >1. In terms of (two-dimensional) 
potential theory language, the singularity of u at Pp is that of a “unit 
point charge.” 

The determination of u and its analytic continuation may be car- 
ried out by reducing it to the case considered in I by the application 
of the conformal transformation 


(3) Z=X+i¥, s=xt+ iy. 


This straightens out the boundary r=1 by transforming it into the 
line X =0; the boundary condition (2) is preserved, while the unit 
point charge at z =e# is carried into unit point charges at 


(4) Z = Ink = H+ 


Proceeding as in I, §5, we find that each one of the periodic array 
of point charges is “reflected” into a point charge and an exponential 
trail of negative charges. When we change back to the z plane, there 
results a positive image at z=1/h and a distributed charge density 


(5) p(x) = — 2ah*x1, for O< x <h. 


In determining the latter, it is to be kept in mind that charge is pre- 
served under conformal transformations, so that p(x) is found from 
p(x)dx = P(X)dX, where P(X) is the charge density along the X axis 
in the Z plane. From I, equation (23), P(X) = —2ae*+», for x< —h. 
As in I, the restriction I, (9), is pertinent. The total amount of 
charge then vanishes so that u vanishes at infinity like O(|2| —1), 

Proceeding to the three-dimensional case and a spherical boundary 
along which (2) holds, we can no longer utilize conformal transforma- 
tions but proceed as follows. 

In order to determine the function u which is harmonic except near 
Po: (x, y, 2) =(h, 0, 0), (4>1), where u becomes infinite like 1/1» (the 
reciprocal distance from Po) while along the unit sphere r = 1 it satis- 
fies the boundary condition (2), start with the expansion 


(6) u = P, (cos 0) + 1/r0, 


n=0 
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for u, (r=1), and utilize the familiar expansions 


(7) [x P,(cos 6) or P,(cos |, 


according as 7 is less than or greater than h, where @ is the angle with 
the x axis and P, are the Legendre polynomials. If (2) is imposed, 


1 (1 + 2a) 


(8) A, = 
Anti (n a 1) 

The first part of A, leads to a point charge of amount 1/h at P;, the 

spherical inverse point of Po; the latter leads to a distributed charge 

density p(x) along the line from the center to P,, provided p(x) satis- 

fies the moment equations 


(9) (2’)(x!) "da! (1 + 2a) 
x x)"4x =— a= 
The solution of these equations is given by 
(1 + 2a)a* 
10 
(10) p(x) re 


Magnification can be used for treating the case where (2) holds 
along a sphere of non-unit radius, while singularities inside the circle 
or sphere can be reduced to outer singularities by means of inversion. 

The two-dimensional case can also be treated in a similar fashion 
with Fourier series replacing the series of Legendre polynomials. 


3. Boundary conditions along two parallel planes. We consider in 
this section continuations of harmonic functions across two parallel 
planes x =0, x = C, corresponding to boundary conditions of the form 
(1) on each plane, though for simplicity we suppose that (1) applies 
along x =0. Then 
Ou 
Ox 


Ou 
(12) = 0, for) 2= iC. 


a 
Cx 


—au=0, for x=0, 


(11) 


Wé shall investigate the reflections of a point charge singularity lo- 

cated at Po: (h, 0, 0) between the planes x=0, x=C, (0<h<C). 
One way of obtaining the singularities of the analytic continuation 

is to apply the familiar method of successive reflections, attending 
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alternately first to x=0, then to x=C. There result point charges 
along the x axis at 


(13) n=0,+1,+2,--- 


? 


as well as proper continuous charge distributions which change ab- 
ruptly to different analytic functions on crossing the points (13). 
This process of reflection gets rather irksome after several stages. 

We shall first follow these reflections by means of the operational 
treatment employed in I, §8. If So(x—A) is the original charge at 
x=h, reflection across x =C leads to the new eharge So(x—2C+h). 
As shown in I, §8, analytic continuation across x =0 results in multi- 
plication by the operator (p+a)/(p—a), in addition to replacement 
of the argument by its negative. This introduces the charges 


pta 


(14) ole +h) + + 2C — 


Further reflection in x=C replaces x by 2C—x, changes f(p) to 
—f(—>p), and results in the charges 


a+ 


(15) ‘ [So(x — 2C — h) + So(x — 4C + h)]. 


Following this by a reflection in x =0 we would apparently have 


ptaa—p 
(16) p—aatp 


[So(x + 2C + h) + So(x + 4C — h)] 
= — [So(x + 2C + h) + So(x + 4C — h)]. 


This, however, may be shown to be incorrect. The explanation of the 
difficulties thus encountered by the operational treatment will appear 
presently. 

We shall give a Bessel-Fourier integral representation for u. Utiliz- 
ing I, equation (16), we put 


(17) f 4008) [ele + + 


where f(A), g(A) are to be determined. Imposing the conditions (11), 
(12) on the integrand and solving, say for f(x), we obtain 


ethh 


A-a 


A+a 


(18) fQ) = 


| 
| 
| 
| 
A-a 
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Expanding f(A) in powers of e~?4*(A—a)/(A+<a) one obtains 


n=0 A+a n=0 A+a 


and the substitution of (19) in the middle term of (17) results in the 
potential function 


=)" 
> edlzth—2(nt dd 


(20) 


> — a\"! 
+ ( ) ed(z-h-2n€) J X add. 
= J, of p) 


To interpret (20) as the potential of charges, note that from I, 
equation (16), it follows that a charge distribution c(x) over the posi- 
tive x axis gives rise to the potential 


(21) f x <0, 
where 
(22) = 


Thus, if f(A) is given (say by means of (18) or (19)), one obtains for 
c(x) a Laplace-Carson integral equation. The solution of the latter 
can be effected by means of a Bromwich integral; in operational no- 
tation the solution appears in the form 


(23) c(x) = f(p)So(x). 


The explicit charge corresponding to the displayed term in the first 
series in (20) is thus 


p—a\" 


The effect of the exponential is to translate the distribution a distance 
2(n+1)C—h in the direction of positive x; there remains to interpret 


(25) = (1 - suey 
= 
This is done by expanding by the binomial theorem and replacing 
(p+a)—*So(x) by x*-'e-**/(k—1)!, for x>0, and by 0 for x <0. Of 


— 
= 
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course, the first term in the expansion of the right-hand side of (25) 
leads to a unit point charge at x =0. It follows that (24) represents a 
point charge at x =2(n+1)C—A and a continuous charge beyond this 
point of charge density P,(x)e~**, where P,(x) is a. polynomial of 
degree n. 

As regards the function g(A) and the last integral in (17), a similar 
treatment is possible resulting in an interpretation in terms of charges 
along the negative x axis. The latter might be described essentially 
as the positive image of the charges for x >C in the plane x = C. How- 
ever, in place of (22) an integral equation of the form 


(26) g(a) = 


now obtains for the charge density c,(x) along x <0. A Bromwich in- 
tegral will not yield the solution of this equation unless one changes 
the independent variable to —x or properly changes the path of in- 
tegration. 

An operational solution of (26) (similar to the solution (23) of (22)) 
could therefore be given only by using a non-orthodox interpretation 
of the operator g(p). Essentially the same situation appeared in I, 
where, in equation (42), 1/(p—a@)So(x) was interpreted as —e**H(— x) 
rather than e*7H,(x), as is usually done in operational calculus* (see 
also I, footnote, p. 884). This also accounts for the difficulties en- 
countered in connection with (16) and the attempted operational 
treatment of the successive reflections (two distinct (and mutually 
inconsistent) interpretations of the operators 1/(p—a), 1/(p+a) oc- 
curring in (15), (16) have to be employed).f 

A very similar treatment can be given for the case wherein general 
linear conditions of the form (11) apply along both x=0 and x=C. 

We believe that some of the features of the above analytic continu- 
ation are characteristic of general differential systems with two point 
boundary conditions. It is planned to touch upon this phase in the 
future. 


GENERAL ELEctTRrIc Co. 
ScHENECTADY, N. Y. 


* Both interpretations yield the Green’s function of the differential equation 
(d/dx —a)u=0 with discontinuity at x =0, the one relative to the boundary condition 
u(+ ©) =0, the other relative to u(— ©) =0. 

+ Roughly speaking, the operational method fails here because of its catering to 
functions which vanish for sufficiently negative x. It can handle the charges for x >0. 
Properly modified (as in I, §8) it can handle the charges for x <0. It is not adapted, 
however, to treating both sets of charges and exhibiting their relations. 


THE SIMPLE GROUP LF(2, 2%) 


ON GENERATING THE SIMPLE GROUP LAF(2, 2”) 
BY TWO OPERATORS OF PERIODS 
TWO AND THREE* 


ABRAHAM SINKOV 


The purpose of this paper is to consider the question of the num- 
ber of abstractly distinct ways in which it is possible for two opera- 
tors of periods two and three to generate the simple group LF(2, 2”). 
The general procedure to be followed in studying such a problem has 
been outlined by Professor Brahanat and has already been applied 
to LF(2, 2%) by the author.{ Those previous results suggested the 
present generalization. 

Since the 2*” —1 substitutions of period two in LF(2, 2”) are all 
conjugate, § it is sufficient in seeking possible generating operators to 
consider only one of them, say 


1,1 
T = (~). 
0,1 
In the representation of G on 2" +1 letters, T leaves fixed the single 


element «. The largest subgroup within which T is invariant is G™, 
composed of all the substitutions 


wf in the GF[2¥], 


which keep the single element ~ unchanged. The group G(*) is 
abelian, of order 2% and of type (1, 1, 1,---). 

If the operators of period three are divided up into sets of com- 
plete conjugates under G“, then the various members of each set 
satisfy with T the same abstract relations, and it is sufficient to select 
from each such set only one operator to serve as a possible second 
generator with 7. The number of these sets to be considered depends 


* Presented to the Society, February 20, 1937. 

t H. R. Brahana, Pairs of generators for the known simple groups whose orders are 
less than one million, Annals of Mathematics, vol. 31 (1930), pp. 542, 543. 

t Necessary and sufficient conditions for generating certain simple groups by two 
operators of periods two and three, American Journal of Mathematics, vol. 49 (1937), 
pp. 69-72, hereafter referred to as Necessary and sufficient conditions. 

§ The various properties of L F(2, 2”) of which use is made in this paper are given 
in Dickson's Linear Groups with an Exposition of the Galois Field Theory, 1901, chap. 
12. 
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on the form of 2%+1. For, when 2" +1 is a multiple of three, G con- 
tains 2"(2" —1) operators of period three, whereas when 2” +1 is not 
a multiple of three, there are 2"(2" +1) such operators. 

In the former case, there are 2” —1 sets of 2" operators each, and 
we may choose as a representative from each set the operator 


0, 
s.=( ), 2" 


where x is a primitive root of the GF[2*]. In the latter case, there 
are 2" +1 sets of 2" operators each. We may again choose the 2" —1 
operators S, together with the two additional operators 


0 
= 
0, 1 


and = (S2v)?. Since 


0, 1 


is of period three, Sev and T generate a tetrahedral group of order 12. 
The same is true of S.v,; and 7, so that neither of these pairs can be 
used to generate the entire group. Moreover, So7 is of period two, so 
that {So, T} is of order six. In every case then, it is sufficient to con- 
sider only the 2” —2 operators S., (a0). 

Now it may happen, for a particular subscript k, that T and S; 
generate a proper subgroup H of G instead of the entire group. We 


note however that 
S.T = 
1, 1+ 


may not be of period 2, 3, or 4 for any a whatever. Hence, on referring 
to the list of possible subgroups of LF(2, 2%), one sees that H is 
either a linear fractional group L F(2, 2¥'*) or else is of order 2*d. But 
the latter case may be ruled out. For, if a subgroup H of order 2*d 
existed, the commutative subgroup within it of order 2* would in- 
clude the commutator subgroup of H. Hence d, and therefore the 
period of S,7, is at most three,* which is not possible. 

H is therefore a linear fractional group, and the number x* must 
be an element in the GF[2¥’"]. If we set N/r=t, then (x*)**!=1, 
from which k(2‘'—1)=0 modulo (2% —1). The number & is thus re- 


*H. R. Brahana, Certain perfect groups generated by two operators of periods two 
and three, American Journal of Mathematics, vol. 50 (1928), p. 348. 
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quired to be a multiple of (2% —1)/(2‘—1). This condition is also 
sufficient, since every subscript which has, with 2” —1, a greatest 
common divisor of the form (2% —1)/(2‘—1), ¢ being a divisor of N, 
corresponds to an operator S which generates with T a subgroup of 
LF(2, 2*) and hence a proper subgroup of LF(2, 2%). 

We may therefore eliminate from the set of operators S, all those 
whose subscripts satisfy the above requirement. There will then re- 
main, in place of the original 2” —2, exactly 


4 QNipip2 — + (— 1)">> 


where the #; are the distinct prime divisors of N and the jth summa- 
tion is taken with respect to all the possible combinations of these #;, 
jat atime. This number A is a simplification of the expression 


(2 — 2) — — 2) + — 2) — - - - 


and its correctness is established as follows: Let k be of the required 
form. Then if ¢ involves 6 distinct primes, S; will appear C;,, ; times in 
the jth summation. The number of times that it will appear in the 
entire expression is 


(1-1) =0. 


Any subscript not of the form required by the preceding paragraph 
will appear only once in the first term. 

The operators S which now remain will all generate with T the en- 
tire group G. But the pairs are not necessarily abstractly distinct, for 
in reducing their number to A we have taken account of only the 
inner automorphisms of G. It still remains necessary to consider the 
possible outer automorphisms. The group of automorphisms of the 
linear fractional group has been studied by O. Schreier and B. L. van 
der Waerden.* It follows from their results that the only possible 
outer automorphisms of LF(2, 2") are to be found among the auto- 
morphisms of the GF[2” ]. From this result, it can be seen that the 
group of outer automorphisms of LF(2, 2”) is the cyclic group gener- 
ated by the substitution 2’ =z, of period N. 

This substitution is commutative with T and transforms S; into 
S21, provided we consider the subscript 2/ reduced modulo (2% —1). 
As a result, the relations satisfied by T and S; are abstractly identical 
with those satisfied by T and Sx,, where 7 may take on any value 


* Abhandlungen aus dem Mathematischen Seminar der Hamburgischen Uni- 
versitat, vol. 6 (1928), pp. 303-322. 
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whatever. Hence the A operators S may be divided up into sets, in 
such a way that the subscripts of all the operators in any set are ob- 
tainable from some one of them by continued multiplications by two 
and reduction modulo (2” —1). The number of operators in any par- 
ticular set is exactly N, except in those cases when the number 2” —1 
and the subscripts of the set have a greatest common divisor of the 
form (2% —1)/(2*—1). But these have already been eliminated. The 
number of sets obtained, and therefore the number of distinct defini- 
tions of the group, is consequently A/N. 

As a corollary, we have the number theoretic result that A/N is al- 
ways an integer. When N is prime, the corollary reduces to a special 
case of Fermat’s theorem. It was pointed out to me by Professor Gill 
that the number A /N is identically Dickson’s Nw,2,* or the number of 
irreducible polynomials of degree N in a GF[2]. This coincidence can 
be used to advantage in the following way: 

Consider two distinct irreducible polynomials in the GF[2] of de- 
gree N in the variables x and X, respectively. Let S, and o; designate 


the substitutions 
0, =) 0, X 
C= 


respectively, and suppose that the pairs of generators T, S; and T, a; 
are abstractly identical. Then we get an automorphism of G by mak- 
ing T correspond to itself and S, correspond to a;. As a result, every 
combination of S; and T will have the same period as the correspond- 
ing combination of a; and T. Consider in particular all the combina- 
tions of the generators which will reduce to the identity. The formal 
procedure of calculating them is identical for the two separate cases 
if we hold all processes of modular reduction by means of the poly- 
nomial modulus until the very last step. We see then that we can thus 
obtain a considerable number of pairs of identical expressions in x 
and X such that both members of each pair will reduce to the same 
value. But such a result will necessarily imply that the two polyno- 
mial moduli are identical, which is a contradiction. Hence a given S. 
will yield with T all the possible definitions of LF(2, 2") if we take 
advantage of all the distinct methods of generating the GF[2¥ ]. 

We shall bring this paper to a close by demonstrating an interest- 
ing general property of every pair of generators S and T of an 
LF(2, 2%). 

In the study of groups generated by two operators of periods two 


* Linear Groups, p. 18. 
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and three respectively, the following substitution has been found very 
useful: 

P = (ST), S = 

Q = (ST)’S, T = PQ. 
It replaces the relations S*=7?=(ST)*=(S-'T-!ST)?=1 by P* 
= Q?=(QP*)?=(QP*)?=1. We shall designate either of these sets of 
relations by the abbreviated notation (2, 3, ; p). 


From the substitution of the equivalents for S and T in the defini- 
tions of P and Q, it results that 


0 (= ) 
By induction 


(= xia ) 
1, x*+ Ay/’ 


where Ay =) ,2-}1x(2*+¥«, Now the period of Q is the same as that of 
the commutator of S and T. Since 


1 Sa 2a 
K = TST-1S-) = 


and 


= (= =) 
1, 1+ By/’ 


where By =1+) it follows that 


and 


1+ By, 1+ By(1 + x**) 


If the period p of K divides 2¥+1, By=x**; if it divides 2% —1, 
By 
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Suppose that p divides 2" +1. Then, since By_, = By +x", we have 
By_,=0, and 


Ay-1 


2*(By_-1 + 1) + = x* + 


2 3 
= 


+1 — ], 
1 


= 1, 
= 1, 


In the same way, if p divides 2" —1, 


(= ) 
= 
1, 4+ 


( 
oF P = 
1 


is of period two. Once again, then, (Q°°+»/?P)?=1. Hence, regardless 
of what particular value p may have, the relation (Q°&+»/?P)?=1 is 
always satisfied. This result has some interesting consequences. 

The relations (2, 3, ; p), (Q°°+»/2P)?=1 define the group desig- 
nated by H. S. M. Coxeter* as G*-*-?. (In this definition, the rela- 
tion Q?=1 is redundantf and may be omitted.) Because of the sym- 
metry of G*-"-, it follows that if n+p, there must exist a pair of 
generators satisfying (2, 3, p; 2), (Q‘"*+»/?P)?=1. Hence, if we con- 
sider all the abstractly distinct definitions of LF(2, 2%) in terms of 
two generators of periods two and three, it follows that the totality 
of values assumed by 2 is identical with the totality of values as- 
sumed by 

Moreover, suppose that by the adjunction of suitable conditions 
we could get a complete definition of LF(2, 2”) in terms of (2, 3, m; p), 
(Q+)/2P)2=1; and suppose further that »¥p. Let us introduce 
the substitution P?=(Q, Q = P? with the additional requirement that 
Pe =1. Then P=Q("t+)/2, and the relations 


* The abstract groups G™-"-?, to appear in the Transactions of this Society. 
t Necessary and sufficient conditions, p. 69. 
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P* = Q? (QP?)? = (QP?) = = (Q(rts)/2pr)2 = 4 


become 
Q” = Pr = (ints) (OP?)3 = (OP)? =1, 


in which the roles of P and Q have been interchanged. We are thus 
enabled to pass directly to a second complete definition in terms of 
(2, 3, p; n), 

It is interesting to observe the first few special cases of LF(2, 2") 
to see how the results of this paper apply. 

When WN =2, the group is the icosahedral group G* *5 which is com- 
pletely defined by the relations (2, 3, 5). This is the only definition; 
hence must be equal to 7, a known result. 

When N=3, the group* is the simple group of order 504. It has 
two definitions, based on (2, 3, 7; 9) and (2, 3, 9; 7). Since Gs =G*-7 9 
it is sufficient in each case to add the one relation (Q&+»/?P)?=1 in 
order to obtain a complete definition of the group. 

No complete definition in terms of two generators has yet been ob- 
tained for any N>3. However, Todd has given a complete defini- 
tiont of LF(2, 2”) in terms of N+2 generators. Since U and R in his 
definition satisfy (2, 3, 2” —1), it follows from the foregoing that these 
two operators suffice to generate the entire group, with the single ex- 
ception of the case N =2, and that one of the values assumed by 7 is 
2" —1. 

When N =4, there are three definitions, and the paired values of n 
and 15, 17; 17: 

As a last instance we mention N =5. Here the six paired values of 
n and pare 31, 31; 31, 31; 33, 33; 33, 33; 31, 11; 11, 31. The duplica- 
tion of the pair 31, 31 means that LF(2, 25) has two distinct abstract 
definitions based on (2, 3, 31; 31). Or, to express it differently, 
LF(2, 25) is obtainable in two abstractly distinct ways as a quotient 
group of G*:31.31. A similar remark holds for the duplication of the 
pair 33, 33. 
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* Necessary and sufficient conditions, p. 70. 
7 Journal of the London Mathematical Society, vol. 11 (1936), p. 106. 
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NOTE ON THE DEFINING RELATIONS FOR THE 
SIMPLE GROUP OF ORDER 660 


F. A. LEWIS 


Dickson* has given the following set of defining relations for the 
simple group of order 660: 


(1) =f, T?=T, (ST)? = TI, 
(2) (S*7S*T)? = I, 
(3) (S‘TS*T)? = I. 
It follows from (3) that 
IST = 
so that 
TST = S*TS'TS"TS'TS§ 
= STS*TS*TS$. 
From this last relation 
SIS IS = I, 
STST = STSTS*TS*, 
(S°TS8T)? = S*TS*TS*TSTS*TS*TS* 
= S*7S*7S' 15° 
m= -S*TS* 
= 
S*7STSTS* 
=f, 
Since (2) has been derived from (1) and (3), it follows that (1) and 
(3) are sufficient to define the simple group of order 660. 
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* This Bulletin, vol. 9 (1903), pp. 204-206. 
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